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SUMMARY 


OBJECTIVE 


Develop an efficient and accurate computational method for rational design of 
frangible nosecaps for air- and surface-launched undersea weapons. 


RESULTS 


The powerful geometry anc finite element model (FEM) pre- and post- 
processor PATRAN1, a potential-flow computer code that can calculate dynamic pres- 
sure-time histories of an arbitrary entry body called ENTRY, and the nonlinear finite 
element analysis (FEA) code ABAQUS have been effectively linked. This linkage 
allows rapid and accurate assessment of the state of stress and deformation of missile 
nosecaps intended to break up at water entry. This technique has been dubbed the 
Water Entry Structural Technique (WEST). 


Two sample nosecap analyses are included in this report. One, a 90-degree 
cone, is an example of a known shape, previously experimentally characterized for 
water entry behavior. This shape serves as the control shape for results evaluation. 
The other is an engineering example of a nosecap that can be used on missiles like 
the vertical launch ASROC (VLA). 


Since no single computer at the Naval Ocean Systems Center (NOSC) has the 
correct combination of processing power and graphics capability to host the various 
parts of WEST, portions of the codes are hosted on various computers at different 
locations within NOSC. The PATRAN pre- and post-processor is hosted on the NOSC 
Code 936-owned VAX 11/785 computer, named FLIPPER, as is the water entry por- 
tion of this code called ENTRY. The nonlinear FEA code ABAQUS is hosted on both 
FLIPPER and the General Purpose Computer Center (GPCC)-owned Convex mini- 
supercomputer STINGRAY. A small computer code called YADAP, used for plotting 
results of pressure-time histories, is hosted on a personal computer (PC). 


An abbreviated operating manual for WEST has been written and is included 
as appendix G. 


CONCLUSIONS 


The objective of this in-house Independent Exploratory Development (IED) 
effort has been met. The code linkage has been validated through comparison with 
experimental work. WEST is a valuable analytical tool that reduces the design cycle 
time for frangible nosecaps used for water entry. 


'PATRAN is a registered trademark of PDA Engincering, Inc. 
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RECOMMENDATIONS 


The code linkage process described in this report is suitable for a scientific 
engineering application within a Naval laboratory such as NOSC. It should not be 
considered a “production” code suitable for commercial application, as there are sev- 
eral “rough edges” within the execution of WEST. In particular, the limitations of the 
water entry portion of the code must be recognized during the development of the 
initial geometry and finite element model (FEM) within the pre-processor portion of 
WEST. While arbitrary shapes can be accommodated in the water entry process, 
there are limitations pertaining to mesh density, spacing, and numbering that pre- 
clude executing any old FEM that happens to describe a particular geometry of a 
nosecap. 


The user must have a background in the PATRAN pre- and post-processor. 
This IED project does not attempt to reduce or foreshorten the complex task of trans- 
lating nosecap drawings into “PATRANese”. A FEA background is essential to create 
viable FEMs that will not only execute throughout WEST, but will truly produce 
“right” answers at the conclusion of the process. As with any FEA technique, it is 
often as much art as science, and not something that can be learned overnight by 
reading this report. 
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INTRODUCTION 


Improvements to the analytical procedure described for the VLA nosecap can 
be made in the following areas: 


(a) ABAPAT translator revision. The PDA Engineering Release 3.0 
ABAPAT translator has numerous “bugs” in it. In particular, while the 
outer fiber normal stresses are computed by ABAQUS and are contained 
in the FILENAME.FIL postprocessing file, ABAPAT cannot translate all 
of them into the FILENAME.STPilj.NOD files needed for complete inves- 
tigation of the results. According to PDA Engineering, Inc. this deficiency 
will be corrected in Version 4.0, which will soon be released. 


(b) Non-linear material model. ABAQUS is equipped to accommodate 
nonlinear stress-strain curves. As the plastic materials typically used in 
nosecaps are nonmetallic, increased accuracy would result by taking 
advantage of this feature. For the work accomplished for this IED project, 
linear-elastic materials properties were used. 


(c) Nonlinear geometry model. In addition to material nonlinearities, 
ABAQUS can accommodate geometric nonlinearities. These should be 
used in the area of the fin-shell interface of the VLA nosecap, to improve 
the accuracy of the relative deformations of the structural elements. 


BACKGROUND 


Early investigations into the phenomena of water entry were concerned with 
the behavior of seaplane hulls upon landing. At the heart of these investigations was 
the determination of the surface pressures, forces, and moments on the entering 
object. Because of the mathematical difficulty of the problem, these investigations 
relied on many assumptions. Rather than attempt to rewrite a compilation of 
approaches found in the literature, an excerpt from Wardlaw, Morrison, and Baldwin 


(1977) follows. 


(NOTE: Superscripted numbers in the following quoted passage refer to the original 
authors’ work. For complete citations see Wardlaw et al. (1977), and the bibliography 
at the end of this report.) 


Attempts to analyze the water-entry problem originate (circa 1929) 
with the work of Von Karman.! Comprehensive surveys of this field are 
provided by May,2 Thingpen,3 Szebehely,4 and Moran.5 The main 
thrust of early work follows the formulation developed by Von Karman 
and Wagner.6 In this approach, a potential flow model is used and 
forces are calculated by the added mass concept. The submerged por- 
tion of the body is often fitted or replaced by another with the same 
surface cross-sectional area for which a closed form solution is avail- 
able. A linearized version of the free-surface boundary conditions is 
applied to determine the surface shape. Most of the theories are 


restricted to vertical entry of simple geometries. In recent years, com- 
putational efforts have been made to obtain a solution using the non- 
linear boundary conditions. An early example of such work is that of 
Chu and Falconer.7 A relaxation method was used to solve the poten- 
tial problem for arbitrary bodies. This project was abandoned due to 
problems with excessive computational time and surface contact dis- 
continuities. The same formulation for the vertical entry of cones has 
been treated by Weber8 using a distribution of source dipoles. More 
recently, Shere and Vander Vorst9 and Vander Vorst and Rogers10 
have used the marker and ceil method to develop a detailed viscous 
model of vertical cone entry. 


Although the description of the pressures, forces, and moments on entry bod- 
ies is of major importance to designers of water-entry bodies, the ability to character- 
ize these loads is only part of the solution to the general engineering design problem 
of frangible nosecaps. Once the water-entry loads are known, the structural perform- 
ance of the entry body must be obtained, and the optimization of the body made to 
accommodate these loads. 


The application of traditional analytical tools to this type of engineering 
design has been supplemented in the last 15 years by the FEA approach. FEA codes 
are currently available, such as MSC/PAL-II,2 CASA/GIFTS,3 and NISA II-PC4 
which run on PCs on a design engineer’s desktop. More powerful versions of these 
codes, such as MSC/NASTRANS execute on higher power workstations, mainframe, 
and supercomputer platforms. 


These codes, however, (especially those hosted on PC platforms) are normally 
used for linear analyses, which are necessarily limited to small displacements and a 
constant relationship between stress and strain in the material. Unfortunately, struc- 
tures designed to fail (frange) upon water impact, such as torpedo nosecaps used in 
weapon systems for air- or surface-launched antisubmarine warfare, cannot be ana- 
lyzed by the linear-elastic method. Due to the nature of the design, these structures 
undergo relatively large strains during the entry phenomena, well beyond the applica- 
tion of linear stress analysis. In addition, these structures are necessarily fabricated 
from some sort of engineering plastic, which typically has an extremely nonlinear 
stress-strain curve. Also, resolution of stresses within complicated nosecap geome- 
tries can result in large problem sizes. Lastly, these structures usually possess some 
form of nonlinear contact boundary condition, which cannot be handled by linear- 
elastic theory. These limitations become insurmountable when dealing with linear 
FEA codes. 


In the early 1980s, NOSC was tasked to design a frangible nosecap for the 
VLA. Because of the highly streamlined profile of the nosecap required for Mach 2 
flight, the subsequent pressure profile generated over the nosecap surface upon water 
entry was significantly smaller than on the antisubmarine rocket (ASROC) nosecap, 


2MSC/PAL-II is a registered trademark of the MacNeal-Schwendler Corporation. 
3CASA/GIFTS is a registered trademark of Computer-Aided Structural Analysis, Inc. 
4NISA II-PC is a registered trademark of Engineering Mechanics Research Corporation. 
SMSC/NASTRAN is a registered trademark of the MacNeal-Schwendler Corporation. 


upon which the initial design was based (Jung, 1984). The VLA nosecap showed a 
marked propensity to remain intact upon water entry, rendering the torpedo payload 


useless. 


A new design concept was subsequently developed by NOSC, in collaboration 
with Loral Systems Group, the prime contractor for VLA. In this design, an ogival- 
shaped plastic shell is filled with segmented rigid foam. The tip of the shell is cut off 
at approximately the 5-inch-diameter station to form a “hollow point” nosecap shell. 
A multipound, lead-filled, aluminum nosetip is attached to the front of the shell. The 
nosetip is designed to remain in place from launch through airframe separation, to 
provide an aerodynamic profile. At parachute deployment, the nosetip is to fall away, 
revealing a high-drag, hollow-nosecap shell designed to break apart at water entry. 
Although complicated by a difficult and restrictive design envelope, a workable design 
was developed. 


In 1988, an effort to find alternative designs was initiated by Jung and Plapp 
(1988). Two designs were investigated: an ogival shell stiffened with tapered foam 
columnar supports, and an ogival shell with blade stiffening. Wardlaw’s, et al. (1977) 
computer code ENTRY was used to determine tne entry body pressure profile. The 
linear analysis capabilities of MSC/NASTRAN were used to generate the subsequent 
stress and displacement time histories of the nosecaps. The analysis revealed two 
shortcomings: First, the manual coding of ENTRY to produce a suitable pressure- 
time distribution on the VLA profile for FEA was extremely laborious and error- 
prone. Second, the linear capabilities of MSC/NASTRAN did not accurately predict 
the structural response of the nosecaps to the dynamic water entry load. Significant 
concerns were voiced by Jung and Plapp (1988) to the effect that potentially viable 
alternatives to the “hollow point” VLA nosecap were being discarded due to the lack 
of suitable analytical tools to properly develop the design of these alternatives. 


This in-house IED project was developed to assuage that concern, and signifi- 
cantly improve NOSC’s capabilities to design frangible nosecaps for water entry. 


APPROACH 


As extensively documented in the literature, the mathematical approach to 
the determination of the water-entry loads on a nosecap is extremely complex com- 
putationally. Therefore, some type of computer code is necessary to solve the frangi- 
ble nosecap design problem. A literature search revealed only a few codes, such as 
PHOENICS6 and PISCES-3DE,7 that can solve fast transient compressible flow phe- 
nomena to characterize the pressure-time history during water entry. Unfortunately, 
these are not public-domain softwares, and are consequently extremely expensive to 
own or lease. As a goal of this project was to obtain an inexpensive solution to the 
water entry structural problem, these codes were not pursued. 


In the late 1970s, an engineering code was developed by the Naval Surface 
Weapons Center/White Oak Laboratory (NSWC/WOL) for calculating pressures and 


SPHOENICS is a registered trademark of CHAM of North America, Inc. 
7PISCES-3DE is a registered trademark of Physics International Company, a California company. 


loads at high speed water entries. Rather than attempt to rewrite in original form the 
authors’ intent, an excerpt from Wardlaw et al. (1977) is included: 


... The present approach differs from such efforts (Von Karman, May, 
Moran, et al.) in several important aspects. Application of Hess and 
Smith’s numerical solution method allows arbitrary bodies to be 
treated and makes it possible to calculate pressures on the surface of 
the body. These pressures are themselves of interest and can be inte- 
grated to provide both force and moment information... ” 


The flow field about the entry body is described by Laplace’s equation. The 
boundary condition 


-Vd:en=VE' en (1) 


is applied on the entry body and ¢ = 0 to the effective planar surface. This surface is 
assumed to rise at the rate (C,- 1) Ve sin (9). Pressures are calculated from succes- 
sive solutions at differing depths using the unsteady Bernoulli’s equation, which is 
cast in a reference frame moving with the model 


Cp = E -ZVe: Vo- cor Ve? (2) 


The preceding equation produces pressure and force coefficients that are 
independent of model scale and entry velocity. The value of these two parameters 
must be simulated through an appropriate choice of the rate of surface rise (i.e., pick- 
ing the correct Cy value). 


To implement the potential flow solution, the nose of the entry body is 
approximated with a series of planar, quadrilateral elements. Each of these elements 
is defined by four points or nodes lying on the body surface. The computation pro- 
ceeds by inserting the model into the water in a series of steps, each at a depth, Ah, 
greater than the previous one. At every step, the group of elements comprising the 
submerged portion of the model is redefined and arranged into a form amenable to 
the potential flow calculative procedure. The nodes defining a particular element are 
checked to determine whether they are above or below the water line. Elements with 
all four nodes above the water surface are discarded, whereas those elements with all 
four nodes below it are included without change. Elements which are intersected by 
the water surface are redefined with an upper edge coincident with it. A constant 
source strength is assumed distributed over the surface of each element, and the 
value of this strength is determined by satisfying eq. 1 at the centroid of each ele- 
ment. The boundary condition, ¢ = 0, is satisfied by locating an image of each ele- 
ment of opposite source strength above the water surface. 


At each depth, the pressure coefficient, C,, is evaluated at element centroids 
using eq. 2, which is a hodv-fixed frame of reference. The quantity V@ appearing in 
this equation is determined directly at each depth, but must be calculated using the 
value of ¢ at the same body location in adjacent steps. In most cases, this require- 
ment does not pose a problem, since ¢ is calculated at the element centroids that are 


at fixed locations on the entry body. Detining ¢-, to be the value of @ at a specific ele- 
ment centroid at step n and using the central differences, 


ag adh ag : Cas ~@ 
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at at dh Vessih cx 6 2AA . | (3) 


The foregoing procedure cannot be applied directly to elements that are intersected 
by the water surface. Under these conditions, the element centroid location on the 
entry body changes slightly from step to step and certain approximates must be intro- 
duced. 


The computational method has been applied to the oblique entry of disk cylin- 
ders, spheres, and to the vertical entry of cones and spheres. In each case C,, is 
assigned a value consistent with existing experimental data. The computed pressure 
coefficients are in good agreement with measurements on the wide variety of bodies 
studied. The calculated pressures reflect not only basic body geometry, but also posi- 
tion on the model surface. The peak pressures tend to be somewhat larger than the 
measured ones. Fortunately, these act over very small areas and thus do not affect 
the calculated loads significantly... 


This computational method is described in detail by Wardlaw et al. (1977), in 
which the computer code ENTRY is described, written in FORTRAN for a Control 
Data Corporation (CDC) computer. As this is public domain software, and thus avail- 
able at no cost to NOSC, it was obtained from NSWC/WOL and installed by Wardlaw 
and Loungs from NSWC/WOL on the NOSC compute: FLIPPER. 


Describing the geometry of the entry body in a form that ENTRY can read can 
be a time-consuming and error-prone task. No graphics pre-processor exists for 
ENTRY. There is no post-processor to enable graphical display of the resulting pres- 
sure-time slices on the geometry either. As it was the intent of this IED task to 
develop a general purpose computational capabiliiy, it was decided iv alter ENTRY to 
enable it to interface with the commercial pre- and post-processor PATRAN. Unfortu- 
nately, PATRAN is not available as public domain software. This compromise was 
accepted upon receipt of assurances from the GPCC for long-term support for 
PATRAN at NOSC. 


There are a number of large, nonlinear FEA codes on the open market such as 
DYCAST/GC,8 and MARC.9 Unfortunately, most of them are extremely expensive to 
obtain, either on a purchase or a lease basis. Some, such as DYNA-3D, are available 
as public domain software, but are so user “unfriendly” as to make them virtually 
useless to the less-than-dedicated user. 


ABAQUS,10 a general purpose FEA program, has special emphasis on non- 
linear structural engineering applications. It has extensive material, element, and 
procedure libraries. Material models include plasticity for metals as well as rubber 


8DYCAST/GC us a regustered trademark of the Grumman Corporate Research Center. 
9MARC is a registered trademark of MARC Analysis Research Corporation. 
1OABAQUS is a registered trademark of Hibbitt, Karlson, & Sorenson, Inc. 


elasticity. It features both user-specified and automatic control of step size and a sim- 
ple, compact problem-definition language. In addition, the David Taylor Research 
Center (DTRC) makes this code available to all Naval laboratories through the Naval 
Engineering Software Support (NESS) Office, for a nominal annual fee. This code 
was chosen for implementation in this IED program. 


Figure 1 describes the linkages and program flow to interconnect PATRAN 
with ENTRY and ABAQUS. Figure 1 also describes the various translators that were 
written to facilitate this linkage. The ability to execute this data flow is the final 
product of this IED project. 


Figure 1 illustrates the flow of information, from the geometry and FEM pre- 
processor PATRAN (black figures) through the water entry code ENTRY (red fig- 
ures), and into the FEA code ABAQUS for stress and displacement analysis (green 
figures), then back again to PATRAN to view the resultant stress and displacement. 


FRANGIBLE NOSECAP ANALYSIS 


To illustrate this process, a step-by-step example for a 90-degree (included 
angle) cone-shaped nosecap is presented (at times referred to by bullet numbers). 


The analysis process begins with the construction of the geometry of the 
nosecap via the PATRAN Phase I modeling technique on FLIPPER. This geometry 
model is then meshed with nodes, and four-noded quadrilateral shell elements to pro- 
duce the Phase II FEM within PATRAN (bullet 1). Figure 2 is a hardcopy image of 
the cone, modeled as a half-symmetric nosecap, created within PATRAN. This image 
shows 100 nodes, and 84 four-noded quadrilateral elements which form the Phase II 
finite element PATRAN model. Figure 3 is a side view of the nosecap. showing the 
90-degree included angle. Figure 4 is an expanded view of the tip of the nosecap, 
showing the four quadrilateral elements used to mesh the flat tip of the nosecap, and 
node 100 at the center of the truncated tip. 


Note that one of the peculiarities of the ENTRY code is that it will not exe- 
cute a pressure-time history on a pointed (zero radius) entry body. Thus, it is neces- 
sary to create a very small flat tip on any pointed geometry to accommodate this 
quirk. 


Upon completion of this Phase II FEM, a request for neutral file generation is 
made within PATRAN. This neutral file contains only the Phase II FEM, and is 
named PATRAN.OUT (bullet 2). This completes the mesh definition portion of the 
problem, and PATRAN is exited. Exiting of PATRAN results in the generation of the 
default PATRAN data file PATRAN.DAT (bullet 3), which must be manually renamed 
to FILENAME.PAT (bullet 4). “FILENAME?” is a user-optional identifier for the par- 
ticular nosecap.) This completes the initial PATRAN pre-processing portion of the 
analysis task. All files generated to this point now reside in the user’s space on 
FLIPPER. 
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Flow chart of frangible nosecap analysis at water entry. 
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Figure 3. Side view of 90-degree blunted cone model. 
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Figure 4. Enlarged view of the tip area of the 90-degree blunted cone. 


Bullet 5 is a translator program written for this project called PATENTR. 
PATENTER resides on, and is executed from the operating system of FLIPPER to 
translate the nodal definition and connectivity generated within PATRAN in neutral 
file format (PATRAN.OUT, bullet 2) into a nodal definition and connectivity format 
readable by the ENTRY code. The output of PATENTR is a file called FILE- 
NAME.LST (bullet 6). Using the FLIPPER text editor, a file containing ENTRY con- 
trol parameters (entry velocity, angle, wetting factor, etc.) is then appended to 
FILENAME.LST and renamed to FILENAME.IN (bullet 7). The sample ENTRY 
input file CONE90.IN is included in appendix A for reference. At this point, ENTRY 
is ready for execution. 


ENTRY is executed from the FLIPPER operating system (bullet 8), using 
FILENAME.IN for the ENTRY input file. This process automatically produces the 
ENTRY default output file FILENAME.OUT (bullet 9), and the file FILENAME.PRE 
(bullet 12), which are saved to the user’s space on FLIPPER. The output file FILE- 
NAME.OUT may be edited by the FLIPPER text editor or printed at the user’s dis- 
cretion. It contains a complete listing of the water entry parameters, the nodal 
definition and connectivity of the geometry, pressure-time history of each element, 
integrated forces, and moments over the entire nosecap. The file CONE90.OUT is 
included in appendix B as an example of a typical ENTRY output file. 


Within ENTRY is an option to produce only the file FILENAME.OUT, and 
exit the entry process (bullet 10). This option would be chosen early in the analysis 
process, for error checking, steady-state drag coefficient calculations, or use of 
ENTRY for some other purpose. If this option is not exercised, then ENTRY uses a 
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portion of the ENTRY code (bullet 11) to produce the file FILENAME.PRE (bullet 
12). This option discards all extraneous information printed in FILENAME.OUT, and 
rearranges the pressure-time history on each element into a format that can be read 
by PATRAN. FILENAME.PRE is created in the user’s space on FLIPPER. 


At bullet 13, a decision is made whether or not to view the pressure distribu- 
tion displayed at the centroids of the elements of the nosecap as a function of time 
during the entry event. If desired, PATRAN is reentered at bullet 14. This execution 
of PATRAN calls upon the PATRAN data file FILENAME.PAT (bullet 4), and the 
formatted results file FILENAME.PRE (bullet 12) to produce color displays of the 
pressure distribution overlaid on the FEM. These graphics ima,es, one per time-step, 
are displayed on a suitable color graphics terminal (bullet 15). These images may be 
captured on hardcopy (bullet 16) if desired. 


Figures 5 through 15 are typical hardcopy images which illustrate this step. 
These figures show the pressure distribution on the 90-degree nosecap during the 
water entry event. Figure 5, at 0.16 ms after vertical entry, shows elements 1 through 
4 fully wetted. The pressure, displayed at the centroid of the quadrilateral elements 1 
through 4, is 24.1 psi. Figure 6, at 0.80 ms after entry, shows the pressure on the flat 
tip rising to 129.0 psi, and the pressure on the ring of elements 5 through 12 wetted 
with a pressure of 131.0 psi at their centroids. Figure 7, 1.61 ms after entry, shows 
the nosecap wetted to the upper edges of elements 13 through 20. The tip pressure 
has decayed to 106.0 psi, while the pressure on elements 5 through 12 has decayed to 
95.7 psi, and the pressure on elements 13 through 20 has risen to 135.0 psi. By fol- 
lowing figures 5 through 15, the progression of the nosecap entry into the water, as 
well as the progression of the pressure distribution up the nosecap, can be seen. At 
figure 15, full submergence has occurred after 8.06 ms. The pressure has decayed 
from a high of 129.0 psi at the tip, to a value of 79.2 psi around the ring of elements 
53 through 60, and has not yet decayed to steady state along the upper ring of 
elements 77 through 84, which have just been wetted. The pressure peak has passed 
up the nosecap, while the earlier-wetted elements have reached their steady state 
drag phase pressure profiles. 


Referring back to figure 1, bullet 17 is a decision whether or not to proceed 
with the FEA portion of this analysis task. If the color hardcopies of the pressure dis- 
tribution, such as figures 5 through 15, are satisfactory to complete the task at hand, 
and no FEA work is intended on the nosecap FEM, then PATRAN may be exited at 
bullet 17 and the task completed. However, if the pressure distributions are accept- 
able to the analyst, and additional FEA work is required, then PATRAN is exited at 
bullet 17 and the WEST process rejoined. 


The translator ENTPRES (bullet 18), another special-purpose translator writ- 
ten for this project, resides on FLIPPER. Execution of ENTPRES from the operating 
system of FLIPPER produces a file called TAPE9.DAT (bullet 19), and the file 
PRES.ABA (bullet 24), both stored in the user’s space on FLIPPER. 
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Figure 10. Pressure distribution on 90-degree blunted cone at t = 4.0298 ms (step 6). 


Figure 12. Pressure distribution on 90-degree blunted cone at t = 5.6417 ms (step 8). 


Figure 14. Pressure distribution on 90-degree blunted cone at t = 7.2536 ms (step 10). 


Figure 15. Pressure distribution on 90-degree blunted cone at t = 8.0595 ms 
(step 11, fully submerged). 


At bullet 20, a decision is made whether or not to view the pressure-time his- 
tory on an elemental basis. If the decision is made to do so, then the user must log off 
FLIPPER and go to the personal computer that hosts YADAP (bullet 21). YADAP is 
an interactive plotting package that plots the pressure on each element as a function 
of time. Through the General Communications Backbone (GCB), a remote login to 
FLIPPER must be made. The file TAPE9.DAT (bullet 9) must then be transmitted 
via the MCP communications protocol from the user’s space on FLIPPER to the 
user’s hard disk on his PC. After executing YADAP, the elemental time-history plots 
can be captured as hardcopy images on a suitable graphics printer or plotter (bullet 
22). 


Figure 16 illustrates this portion of the analytical process, showing pressure- 
time histories for “rings” or “groups” of elements. Element group 1 consists of ele- 
ments | through 4 on the flat face, element group 2 is the ring of elements 5 through 
12, and so forth. These YADAP plots show, on an elemental basis, how the pressure 
peaks at a later time the farther up the nosecap the element is located, as well as the 
reduction in peak pressure the farther from the nosecap tip the element is located. 


Another decision point is reached at bullet 23. Here, if the analyst decides 
that there is enough information about the problem, he/she can exit the process. 
However, if the analyst is satisfied with the pressure distribution time-histories to 
date, (figure 16) but wishes to do the FEA on the nosecap, he/she returns to the main 
program flow in the program stream. 
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Figure 16. Pressure-time histories (0.000-0.008 s). 
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Earlier in the process, the translator ENTPRES (bullet 18) produced a file 
called PRES.ABA (bullet 24), which was stored in the user’s space on FLIPPER. This 
file is a portion of the “history” input file needed by ABAQUS (bullet 25) for the FEA 
of the nosecap. To execute ABAQUS and produce the FEA of the nosecap, it is neces- 
sary to execute the translator PATABA (a PDA Engineering translator, part of the 
PATRAN software) from FLIPPER (bullet 26). This translator operates on the 
PATRAN neutral file PATRAN.OUT (bullet 2) previously created and stored on 
FLIPPER, and produces the ABAQUS model input file called ABAQUS.INP (bullet 
27). This file must then be manually renamed to FILENAME.INP (bullet 28). After 
appending it to the file PRES.ABA (bullet 24), it is possible to complete the FEA por- 
tion of this process by executing ABAQUS in bullet 25. 


ABAQUS is executed from the operating system of FLIPPER, using FILE- 
NAME. INP as the input file. Successful execution of ABAQUS in bullet 25 on FLIP- 
PER produces several permanent files of the analyst’s choice (bullet 29), and a myriad 
of temporary files which are erased at program halt. These files are stored in the us- 
er’s space on FLIPPER. Typically, the following are the files which are retained at 


this point: 


FILENAME.INP input file 
FILENAME.DAT output file 
FILENAME.FIL post-processing output file 
FILENAME.LOG accounting file 


The sample ABAQUS input file for the 90-degree cone (CONE90.INP) is 
shown in appendix C. 


Upon completion of ABAQUS, the translator ABAPAT (another PDA Engi- 
neering translator, bullet 30) is called upon to operate on the post-processing output 
file FILENAME.FIL (bullet 29). ABAPAT produces the PATRAN displacement, ele- 
mental stress, and average nodal stress files (bullet 31) called STPilj.DIS, STPilj. 
ELS, and STPilj.NOD respectively. In these files, i is the step number, and j is the 
increment number within each step. Obviously, a complex nonlinear dynamics prob- 
lem can have many steps, and numerous increments within each step. Thus, this por- 
tion of the code tends to produce massive amounts of data, that the analyst is well 
advised to prescreen for need. 


Once the displacement and stress files are complete in bullet 31, the post- 
processing of the results can begin. In bullet 32, after executing ABAPAT, the 
PATRAN environment is re-entered. Here, after calling up the PATRAN model file 
FILENAME.PAT (bullet 4), displaced geometry and stress contours for each time 
increment during the entry process can be created. These images can be viewed on 
the appropriate graphics terminal (bullet 33) or output as color hardcopy in bullet 34. 


For example, figure 17 is a hardcopy image of the displacement of the 
90-degree nosecap 0.15 ms after entry. As the tip is just wetted, no appreciable defor- 
mation of the nosecap structure has occurred. Figure 18, at 0.80 ms after entry, 
shows the beginnings of deformation of the nosecap under the pressure distribution 
previously illustrated at this same time-step in figure 6. By following the progression 
from figures 17 through 23, the deformation of the nosecap in response to the pres- 
sure-time history can be clearly seen. This deformation is greatly enlarged for clarity. 
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= 1.6119 ms (step 3). 
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Figure 19. Deformed shape of 90-degree blunted cone at t 
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2.4179 ms (step 4). 
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Figure 20. Deformed shape of 90-degree blunted cone at t 
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Figure 23. Deformed shape of 90-degree blunted cone at t = 4.8357 ms (step 7). 


Deformations of elastic bodies result in associated stresses. Figure 24, a maxi- 
mum distortion energy theory (Von Mises) stress plot, shows the state of stress of the 
nosecap at time-step 2, 0.80 ms after water entry. As only the tip of the nosecap has 
been wetted, the Von Mises stresses have only built to a low value on the nosecap. 
Note that, since no effort was made to accurately model the upper supported edge of 
this nosecap mode, the actual stress values reported in the stress spectrum along the 
right edge of the figure are meaningless. They serve only as relative values during the 
event. 


Figures 25 and 26, time-step 3 at 1.61 ms after entry, indicate that the stress 
levels in the blunted tip have risen to a level exceeding 2.72 x 104 psi. 


Figures 27 through 30 are a progression of the Von Mises stresses through the 
nosecap during the water-entry process, at the same time-steps as the previously dis- 
played deformation histories (figures 17-23). The stress levels along the edge can be 
seen to fall during the entry process, and then rise again as the stress wave 
approaches this boundary. As the stress wave follows the pressure pulse up the shell, 
the boundary condition imposes an inordinate amount of bending into the shell at the 
upper edge. This bending skews the stress contouring to the point where the analysis 
is considered unrealistic beyond the point where the nosecap is fully wetted. 


As this nosecap has been included only to illustrate the flow of information 
through this analytical process, the final few figures indicate some of the details of 
this discipline that require extreme attention by the analyst. Blind faith in computer- 
generated numbers easily leads to unrealistic expectations from the FEA approach 
and dangerously incorrect answers. 


Assuming the analyst is satisfied with the progress to date, production of 
these displacement and stress contours completes the program flow at bullet 33. The 
propensity of the nosecap to break up at some point in the entry sequence can be 
determined at this time and a decision made on optimization of the nosecap struc- 
ture. Obviously, if the results are not satisfactory, the program can be started over 
(bullet 35), or exited at this point. 


RESULTS 


Jung and Plapp’s (1988) technical report (TR 1221) dealt with a frangible 
nosecap design for the VLA program. In that document, it was reported that a frangi- 
ble one-piece nosecap was not suitable for the current VLA, due to the difficult, 
restrictive design envelope. 


The approach used to develop this conclusion was a manual design optimiza- 
tion of the nosecap shell. The water entry pressure profile was generated using 
ENTRY, and the resultant stress analysis completed using the linear-elastic capabili- 
ties of the MacNeal/Schwendler version of MSC/NASTRAN. 


To further illustrate the engineering application of the WEST developed in the 
current IED project, it was decided to review the water entry portion of the VLA 
nosecap analysis. 


Figure 31 is the nodal pattern of a half-symmetric 1.25-caliber Von Karman 
ogival nosecap water entry FEM, suitable for use on the VLA missile. This nosecap 
consists of 262 nodes on the exterior of the shell. Figure 32 is a plot of the 228 four- 
noded quadrilateral elements which form the elemental connectivity for the water 
entry model, while figure 33 shows an enlarged view of the tip of the nosecap FEM. 
These figures were produced within the PATRAN portion of WEST (bullet 1, fig- 
ure 1). 


Appendix D is the resultant ENTRY input, entitled CAP90.IN;6. This is the 
result of bullet 7 in figure 1 and is used in the execution of the ENTRY portion of 
WEST in bullet 8. 


Appendix E is the hardcopy output from ENTRY (bullet 9a). Figures 34 
through 40 show the plotted resultant pressure distribution as a function of time for 
the ogival nosecap undergoing vertical entry at a velocity of 130.0 ft/s during the 
time period 0.0 to 3.92 ms. Note that only even time-steps after step 2 are shown for 
brevity. 


These pressure-time histories are shown as element groups (elements grouped 
as rings in the axial direction) in figure 41. This figure shows how the pressure peak 
decays in time during water entry. 
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Figure 31. Nodai point pattern of half-symmetric 1.25-cal. Von Karman ogival nosecap 
water-entry model. 
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Figure 32. Element-numbering pattern of half-symmetric 1.25-cal. Von Karman ogival 
nosecap water-entry model. 
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Figure 33. Enlarged view of the tip area of ogive model. 
Figure 34. Vertical water-entry nosecap impacting water at 130 ft/s, pressure distribution 


at t = 0.0445 ms (step 1, initial wetting). 


Figure 35. Pressure distribution on vertical water-entry nosecap at t = 0.1770 ms (step 2). 
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Figure 36. Pressure distribution on vertical water-entry nosecap at t 
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1.5624 ms (step 6). 
2.3226 ms (step 8). 
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Pressure distribution on vertical water-entry nosecap at t 


Figure 37. Pressure distribution on vertical water-entry nosecap at t 


Figure 38. 
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Figure 39. Pressure distribution on vertical water-entry nosecap at t 
Figure 40. Pressure distribution on vertical water-entry nosecap at t 
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Figure 41. Pressure-time histories (0.000-0.010 s). 
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To increase computational efficiency, it is possible to take advantage of sym- 
metry in the FEM generated in bullet 1. Figure 42 is a nodal plot of a quarter- 
symmetric FEM of the nosecap. Notice that this FEM includes two “blade stiffeners” 
inside the nosecap, which are also composed of four-noded quadrilateral and three- 
noded shell elements, none of which are wetted during the water entry process. By 
re-entering PATRAN, it is possible to modify the figures 31 and 32 FEMs to reduce 
the nodal and elemental count, and add the internal stiffeners, all without compro- 
mising the pressure distribution information generated to date. As this is a non- 
generic detail of the FEM process, it is not specifically shown in figure 1, but can be 
accommodated by the experienced user in the WEST process. In addition, the elemen- 
tal resolution in the circumferential direction has been increased from four to eight 
elements per 90-degree segment. As reversed bending is expected in the unsupported 
regions in the cen‘er of the nosecap, this increased resolution will improve the com- 
putational accuracy of the FEA. 


The figure 42 FEM, after being transferred to neutral file format in bullet 2, 
translated into ABAQUS input format in bullet 26, and after being appended with the 
pressure-distribution information generated in bullet 24, was used to produce the 
ABAQUS input file CAP90.INP;10. This input deck is included as appendix F of this 


report. 


Execution of ABAQUS in bullet 25 produces the output files CAP90.DAT, 
CAP90.FIL, CAP90.LIS, etc. The post-processing file CAP90.FIL is operated on by 
the results translator, ABAPAT, in bullet 30 to produce the hardcopy deformation 
plots shown in figures 43 through 52, from 0.04 ms to 6.26 ms after water entry. 


Obviously, after the initial wetting shown in figure 43, negligible deformation 
has taken place, despite the extreme pressure on the tip of the nosecap (see figure 34). 


The progression of deformation in the nosecap shell with respect to time dur- 
ing water entry is shown in figures 43 through 52. Notice that a “hump” develops just 
upstream from the two blade stiffeners, which is first visible in figure 46 at 1.56 ms 
after water entry. In addition, a concavity in the shell forms between the stiffeners, 
first visible in figure 47 at 2.32 ms after entry. These deformations become progres- 
sively more severe the further the nosecap is immersed. 


As previously illustrated with the 90-degree cone example, deformation of a 
structure body is accompanied by associated stress levels within the structure. In this 
stress analysis process of the nosecap, the nosecap shell is predicted to fail as the 
stress in the shell material reaches some measure of a failure criterion. In the case of 
nonmetallic brittle materials in a nonuniform geometrical arrangement such as a 
nosecap, it is appropriate to use the maximum principal stress theory for a failure 
criterion. This theory asserts that failure or fracture of a material occurs when the 
maximum principal stress at any point in the structure reaches a critical value 
regardless of the other stresses. The critical value of stress is usually determined in a 
tensile experiment, where the failure of a specimen is defined to be either excessively 
large elongation or fracture. Failure is characterized by the separation, or the cleav- 
age, fracture. This mechanism of failure differs drastically from a ductile fracture 
typical of metallic materials, which is accompanied by large deformations due to slips 
along the planes of maximum shearing stress. 
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Figure 42. Nodai point pattern of a nosecap quarter-symmetric model with internal 
blade supports. 


Figure 43. Progression of deformation of the nosecap shell, at t = 0.0445 ms 
(step 1, initial wetting). 
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Figure 45. Deformed shape of shell at t = 0.8440 ms (step 4). 
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Figure 47. Deformed shape of shell at t = 2.3226 ms (step 8). 


48 


Mg 
gee 


Figure 49. Deformed shape of shell at t = 3.9200 ms (step 12). 
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Figure 51. Deformed shape of shell at t = 5.5041 ms (step 16). 
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Figure 52. Deformed shape of shell at t = 6.2642 ms (step 18). 


The ultimate tensile strength (UTS) of the fiber-reinforced PHENOLIC mate- 
rial for the VLA nosecap has an average UTS in the transverse (hoop) direction of the 
nosecap of 5.68 x 103 psi (5.68 KSI) and a UTS of 7.57 KSI in the longitudinal (axial) 
direction of the nosecap. Therefore, failure of the nosecap can be predicted any time 
the value of any of the three principal stresses exceed these values. Note that, as this 
is a thin shell, the through-the-thickness stress gradient is virtually zero, and will not 
be reported. 


Figures 53 through 62 are stress contour plots of the inner fiber first principal 
stress of the nosecap, from initial entry to 6.26 ms after entry. Nodal values are 
shown in the color bar to the right of the plots, which are averaged by PATRAN 
throughout the elements to produce the color contours shown. To reduce repetition, 
only five time-steps are reported. 


Figure 53, at 0.04 ms after entry, shows the maximum first principal stress 
occurring in the second and third rings of elements from the tip. As the nosecap is 
barely wetted at this time-step (see figure 35), the deformation has not had a chance 
to propagate to any great extent, and the stress level has only risen to -4303 psi. 


Figure 54 indicates that the maximum stress has risen to -4565 psi in several 
areas along the midspan between the tip and the stiffener interface. The concavity in 
the shell along the center meridian has started to develop significant stresses of its 
own. 


Figure 55, at 3.11 ms after entry, indicates that the stress level in the midspan 
area along the edges of the FEM have reached a stress level of -11,648 psi, while the 
center meridional stress level has increased to approximately -5400 psi. 
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In figure 56, at 4.73 ms after entry, the stress in the midspan along the edges 
has fallen to -11,327 psi, with a corresponding decrease to -5286 psi along the center 
meridian. The point of maximum siress las moved up the nosecap in the approximate 
location of the water surface. 


Figure 57, at 6.26 ms after entry, indicates that the midspan stress has contin- 
ued to fall, along with the meridional stress. As previously shown in figure 41, the 
magnitude of the pressure distribution decreases with time as the nosecap enters the 
water. As the shell diameter increases with distance from the tip, the stress level 
would be expected to increase with time, had the pressure distribution been constant 
with time. However, this is obviously not the case, and the stress level in the shell 
has been shown to peak at some time after entry, but before the shell is fully wetted. 


To reduce redundancy, the second principal stress distribution is not reported. 


Figures 58 through 62 are stress contour plots of the third principal stresses, 
from 0.18 to 6.26 ms after water entry. In general, these follow the same trends as 
reported in figures 53 through 57. However, the third principal stresses are all tensile 
(positive in magnitude) stresses, and are greatly influenced by the reverse bending 
that occurs in the equatorial direction around the shell. 


A cross-sectional view of the nosecap along the axial meridians reveals some 
bending behavior occurring near the “hump” in the shell. However, the predominant 
deformation is compressive, yielding predominantly compressive normal stresses in 
the axial direction. A cross-sectional view of the nosecap along the hoop direction 
equators reveals considerable reverse bending response in the hoop direction. It is 
this bending response that dominates the stress field used to compute the third prin- 
cipal stresses. 


As the direction cosines of the normal stresses are not reported by ABAQUS, 
it is extremely difficult to describe the directions of the principal stresses. Although 
care was taken in this analysis to insure that the elemental numbering scheme used 
resulted in a uniformly “outboard normal” for each element, ABAQUS sorts the resul- 
tant principal stresses by magnitude. Therefore, it is not possible to state that the 
first principal stresses are aligned in predominantly axial direction, or that the third 
principal stresses are predominantly aligned in the hoop direction. In addition, the 
deformation in response to the pressure-time history is extremely complex, both due 
to the changing loads, and the presence of the blade stiffeners. The principal stresses 
are simply compared to UTS in the maximum principal stress theory. If a failure is 
indicated in one of the three principal directions, it is not important to know in what 
direction that stress lies. 


Figures 53 through 62 show that a dynamic state of stress exists within the 
nosecap during water entry. To attempt to simplify the large amount of data col- 
lected, figures 63 through 66 were condensed from the data. These are plots of the 
state of stress as a function of time for four points on the nosecap shell, aligned with 
the 90-degree meridian and the blade stiffener. Node 28 is located near the tip of the 
shell, node 73 in the midspan between the tip and the blade-shell interface, node 118 
near the lower edge of the blade stiffener, and node 136 near the middle of the blade- 
shell interface (see figure 42). 
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These stress-time histories indicate an obvious periodicity of the data. This 
could be indicative of the period of the moving stress wave within the shell, or the 
lowest order response of the shell. It is also likely that there are considerable higher 
frequency stress waves propagating through the structure that either the nodal spac- 
ing in the FEM was unable to resolve, or the time steps used in the dynamic analysis 
were too large. In any case, there is no need to lend any large degree of engineering 
significance to this periodicity. 


Figures 63 and 64 are the plots of the first principal stresses at nodes 28, 73, 
118, and 136 on the inner and outer fibers respectively. Note that the inner fiber 
stress reaches a maximum value of -12,560 psi at 4.32 ms after water entry, while the 
outer fiber first principal stress for the same node has a value of 0 psi at the same 
time. This apparent anomaly is due to the method that ABAQUS uses to sort princi- 
pal stresses. Although not shown, it is suspected that the second principal stress at 
this time period and nodal location has a nonzero stress value of some negative value, 
which has a great deal more engineering significance than the zero value shown at 
node 118 at this time-step. 


Figures 65 and 66 show the third principal stresses for the same nodes on the 
inner and outer surfaces respectively. In figure 65, it is seen that the third principal 
stress at both nodes 118 and 136 exceed the value of UTS shortly after 3.0 ms after 
water entry, reaching a peak value of 7312 psi. 


At the same time, the third principal stress on the outer surface in figure 66 
reaches a maximum value of 21,240 psi at node 136, while the value at node 118 rises 
to 12,650 psi. Each of these values greatly exceed the UTS of 5680 psi, shown as a 
horizontal line in figures 65 and 66. 


Figure 67 is a maximum envelope of the stress-time history of the entire 
nosecap. This indicates the maximum value of the second principal stress (17,856 
psi), and the maximum value of the third principal stress (21,129 psi) occur simulta- 
neously at 3.11 ms after water entry. The corresponding minimum values are shown 
in figure 68. 


The UTS of the material is 5.68 KSI (Dimelfi, July 1988), while the ultimate 
compressive strength (UCS) is 30.16 KSI (Fogg, October 1989). Therefore, the mini- 
mum margin of safety of the nosecap during water entry occurs in tension at 3.11 ms 
after entry, and has a value of: 


MS (UTS/Max Principal Stress)-1 
(5680/21,129)-1 
~0.73 


A dynamic stress analysis of this nosecap structure is a very complicated 
proposition. Not only are the stress levels affected by the blade-shell interface, but 
the presence of a dynamic load, with associated moving stress waves and the possibil- 
ity of high-frequency shell responses, makes interpretation of these stress contours 
very difficult indeed. As this nosecap example was included to illustrate the WEST 
process, little attempt was made to exactly model the nosecap structure. Thus, while 
these results are considerably more accurate than those reported by Jung and Plapp 
(1988), they should only be used as an indicator of the ability of the WEST to aid the 
analytical process used to design this type of frangible nosecap. 
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CONCLUSIONS 


The objective of this IED effort has been met. The WEST has been validated 
through two separate examples of nosecaps, a 90-degree cone and an engineering 
example of a typical VLA nosecap. While both these examples described undergoing 
vertical entry, the WEST is equally facile at oblique entry. Note that the use of sym- 
metrical FEM must be adjusted accordingly, as an oblique entry results in nonsym- 
metric loading. In addition, a greatly increased amount of data results from the 
analysis, as each element has a unique pressure-time history throughout the water 
entry process. 


Development of the WEST has led to a many-fold reduction in the time 
required to execute a meaningful design iteration of a frangible nosecap. In addition, 
due to the ability to use both geometric and material nonlinearities, the potential 
accuracy of the solution has greatly improved as well. Note that the analytical tech- 
niques described in this report did not take advantage of these features. 


While investigation of an alternate frangible nosecap for the VLA project is 
certainly not complete, the results from the water entry portion of that task indicate 
that breakup at water entry appears likely. 
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Figure 63. VLA nosecap first principal stresses of the inner surface. 
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Figure 64. VLA nosecap first principal stresses of the outer surface. 
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Figure 65. VLA nosecap third principal stresses of the inner surface. 
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Figure 66. VLA nosecap third principal stresses of the outer surface. 
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Figure 68. VLA nosecap minimum principal stresses. 
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RECOMMENDATIONS 


The code linkage process described in this report is only suitable for a scien- 
tific engineering application within a Naval laboratory such as NOSC. It should not 
be considered a “production” code suitable for commercial application, as there are 
several “rough edges” within the execution of WEST. In particular, care must be 
taken during the development of the initial geometry and FEM within the pre- 
processor portion of WEST, to acknowledge the limitations of the water entry portion 
of the code. While arbitrary shapes can be accommodated in the water entry process, 
there are limitations pertaining to mesh density, spacing, and numbering that pre- 
clude executing any old FEM that happens to describe a particular geometry of a 
nosecap. 


In addition, it is assumed that the user has a background in using the 
PATRAN pre- and post-processor. This IED project does not attempt to reduce or 
foreshorten the complex task of translating nosecap drawings into “PATRANese”. A 
FEA background is essential to create viable FEMs that not only will execute 
throughout WEST, but will truly produce “right” answers at the conclusion of the 
process. As with any FEA technique, it is often as much art as science, and not some- 
thing that can be learned overnight by reading this report. 


Improvements to the analytical procedure described for the VLA nosecap can 
be made in the following areas: 


(a) ABAPAT translator revision. The PDA Engineering Release 3.0 
ABAPAT translator has numerous “bugs” in it. In particular, while the 
outer fiber normal stresses are computed by ABAQUS and are contained 
in the FILENAME.FIL post-processing file, ABAPAT cannot translate all 
of them into the FILENAME.STPilj.NOD files needed for complete inves- 
tigation of the results. Therefore, they are not included in this report. 
According to PDA Engineering, this will be fixed in Release 4.0, which is 
due for release within 90 days. 


(b) Nonlinear material model. ABAQUS is equipped to accommodate non- 
linear materials models. As the plastic materials typically used in 
nosecaps exhibit considerably nonlinear stress-strain behavior, increased 
accuracy would result by taking advantage of this feature. For the work 
accomplished for this IED project, linear-elastic materials properties were 
used. 


(c) Nonlinear geometry model. In addition to material nonlinearities, 
ABAQUS can accommodate geometric non-linearities. These should be 
used in the area of the fin-shell interface of the VLA nosecap, to improve 
the accuracy of the relative deformations of the structural elements. 


(d) Mesh rezoning. As the nosecap deforms, the geometry presented to the 
water surface during water entry changes. Yet, in the present form of the 
WEST process, the pressure-time histories are considered to be developed 
by a constant, invariant geometry. With significant deformation, this can 
introduce significant error in the nrocess. A worthwhile improvement 
would be to halt the deformation } - -cess at some point in the analysis 
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loop, remesh the deformed geometry, use this new meshed geometry as 
input to the ENTRY code, re-apply the revised pressure-time histories to 
the deformed geometry, and re-start the analysis. This iteration could be 
performed several times if necessary, until the analyst was satisfied that 
accuracy had been achieved. 


Investigation of the VLA nosecap shape during the development of this WEST 
revealed several areas for improvement in the analytical process used to examine the 
nosecap structure. Because of the negative margin of safety at water entry, additional 
investigation into the design optimization. of this nosecap should be pursued, followed 
by fabrication and test of sample one-piece nosecaps. The potential for success is 
there, and the rewards for major simplifications to the nosecap design are many. 
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ABAQUS 


ASROC 
ABAPAT 
CASA/GIFTS 


Cy 

Ah 

DTRC 
DYCAST/GC 


DYNA-3D 
ENTPRES 
ENTRY 


Pen 

FEA 

FEM 
FLIPPER 
FORTRAN 
GCB 


GLOSSARY 


An FEA computer code that is a registered trademark of Hibbitt, 
Karlson, & Sorenson, Inc. 


Antisubmarine rocket 
A PDA Engineering translator 


An FEA computer code that is a trademark of Computer-Aided 
Structural Analysis, Inc. 


Control Data Corporation 

Pressure coefficient 

The gradient operator 

Wetting factor 

Increment in effective depth between successive steps 
David Taylor Research Center 


An FEA computer code that is a registered trademark of the 
Grumman Corporate Research Center 


An FEA computer code 
A NOSC-developed translator 


A computer code for the prediction of impact loads of an arbitrary 
body during water entry 


The value of ¢ at a specific element centroid at step n 
Finite element analysis 

Finite element model 

A NOSC Code 936-owned VAX 11/785 computer 

A computer language 

General Communications Backbone 

General Purpose Computer Center 

Independent Exploratory Development (Program) 
1000 pounds per square inch 

A velocity potential function 


An FEA computer code that is a trademark of MARC Analysis 
Research Corporation 
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ms 


MSC/NASTRAN 


MSC/PAL-II 


MCP 

MS 

NESS 
NISAII-PC 


NOSC 
NSWC/WOL 
PC 
PATRAN 
PHOENICS 


PISCES-3DE 


PATABA 
PATENTR 
psi 
STINGRAY 
UCS 


Millisecond 


An FEA computer code that is a registered trademark of the 
MacNeal/Schwendler Corporation 


An FEA code that is a registered trademark of MacNeal/ 
Schwendler Corporation 


A VAX command to download a VAX file to PC, or vice versa. 
Margin of safety 
Naval Engineering Software Support (Office) 


An FEA computer code that is a registered trademark of 
Engineering Mechanics Research Corporation 


Naval Ocean Systems Center 

Naval Surface Weapons Center, White Oak Laboratory 
Personal computer 

A registered trademark of PDA Engineering, Inc. 


An FEA computer code that is a registered trademark of CHAM 
of North America, Inc. 


An FEA computer code that is a trademark of Physics 
International Company 


Another PDA Engineering translator 

A NOSC-developed translator 

Pounds per square inch 

A GPCC-owned convex mini-supercomputer 

Ultimate compressive strength 

Ultimate tensile strength 

Vertical launch ASROC 

Water Entry Structure Technique 

A PC computer code for display and plotting of time-histories 
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0.155475445E+02, Q.45849890 , 0.458498 1 
Q0.155475445E+02, 0.362535143E+01, 0.537676239E+01 
0.155475445E+02, 0.252392459E+O1, 0.597459745E+01 
0.155475445E+02, 0.130678225E+01, Q.635241842E+01 
Q.155475445E+02, Q.QQQQQ0000E+00, 0.648415327E+01 
Q.163800011E+02, Q0.€55000067E+01, O.Q0Q0QQ0000E+00 
Q.163800011F+02, 0.641692781E+01, Q.132005358E+01 
Q.1638000 , ©.603526926E+01, Q.254955578E+01 
Q.163800030 , 9.543136311E+01, 0.366216779E+01 
0.163800030E+02, 0.463154984E+01, 0.463154984E+01 
0.163800030E+02, 0.366216731E+01, 0.543136358E+01 
0.163800030E+02, 0.254955530E+01, Q.603526926E+O1 
Q.163800011F+02, 0.132005262E+01, 0.641692781E+01 
Q.163800011E+02, O.Q0QQQQ0000E+00, O.655000019E+01 
Q.171679993E+02, Q0.655000019E+01, 0.O0Q0000000E+00 
0.171679993E+02, 0.641692734E+01, 0.132005370E+01 
Q.171679993E+02, 0.603526878E+01, 0.254955602E+01 
0.171679993E+02, Q0.543136311E+01, 0.366216755E+0 
0.171679993E+02, 0.463154984 0.463154984E+0 
0.171679993E+02, 0.366216803E+01, 0.543136358E+01 
Q.171679993E+02, 0.25495562 1, Q.603526926E+01 
0.171679993E+02, 0.132005453E+O01, 0.641692877E+0 
0.171679993E+02, Q.QOQQQOQ00E+00, 0.655000114E+01 
0.179559975E+02, Q.655000019E+O1, 0.Q00Q00Q0Q000E+00 
0.179559975E+02, 0.641692734E+01, Q.132005370E+01 
0.179559975E+02, 0.603526878E+01, 0.254955602E+01 
0.179559975E+02, Q0.543136311E+01, 0.366216755E+01 
0.179559975E+02, Q.463154984E+01, 0.463154984E+01 
0.179559975E+02, Q0.366216779E+O1, 0.543136358E+01 
Q.179559975E+02, 0.254955626E+01, 0.603526926E+01 
0.179559975E+02, Q.132005453E+01, 0.641692877E+0 
0.179559975E+02, O.QOQQQQQ00E+00, 0.655000114E+01 
0.187439976E+02, O0.655000019E+O01, 0.00Q000Q000E+00 
0.187439976E+02, 0.641692734E+01, 0.132005370E+01 
Q.187439976E+02, Q0.603526878E+01, Q.254955602E+01 
Q.187439976E+02, 0.543136311F+01, 0.366216755E+01 
0.187439976E+02, Q.463154936E+0 0.463154984E+0 
0.187439976E+02, 0.366216731E+01, 0.543136358E+01 
0.187439976E+02, Q.254955578E+01, 0.603526926E+01 
0.187439976E+02, 0.132005358E+01, 0.641692877E+01 
0.187439976E+02, O.QQQQQQ0Q00E+00, O.655000114E+01 
0.19531997 2, 0.655000019E+01, 0.Q0QQ0Q0Q000E+00 
Q.195319977E+02, 0.641692734E+01, 0.132005370E+01 
0.195319977E+02, Q.603526878E+01, 0.254955602E+01 
Q.195319977E+02, 0.5431363 1, Q.366216755E 
Q.195319977E+02, Q.463154936E+01, 0.463154984 
Q.195319977E+02, 0.36 521673 1, 0.543136358E+01 
0.195319977E+02, 0.254955578E+O1, 0.603526926E+01 
0.195319977E+02, O.132005358E+01, 0.64169287 i 
Q.195319977E+02, Q.QQQQQQQ00E+00, Q.655000114E+O1 
0.203199997E+02, 0.655000019E+01, 0.Q00000Q000E+00 
Q.203199997E+02, 0.641692734E+01, Q.132005370E+01 
0.203199997E+02, 0.603526878E+O1, 0.254955602E+01 
0.203199997E+02, 0.543136311E+01, 0.366216755E+01 
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MERID{ 1, MERID22 /MERID34 , MERTID45 7 
: PLL : ,»MERTD56 ,MERID67 , MERID79 


324, 0.540338612E+01, 
325, 0.608881330E+01, 
326, 0:677424002E+01, 
327, 0.745966768E+01, 
328, 0.814509487E+01, 
329, 0.883052254E+01, 
330, 0.267233920E+01, 
331, 0.335655332E+01, 
332, 0.404076719E+01 
333, 0.472498083E+01, 
334, 0.54091949 ; 
335, 0.60934090 ij 
336, 0.677762222E+01, 
337, 0.74618368 : 
338, 0.814605045E+01, 
339, 0.883026505E+01, 
340, 0.268300080E+01, 
341, 0.336600137E+01; 
342, 0.404900217E+01, 
343, 0.473200274E+01, 
344, 0.541500330E+01, 
345, 0.609800434E+01, 
346, 0.678100443E+01, 
347, 0.746400595E+01, 
348, 0.814700603E+01, 
349, 0.883000755E401, 
350, 0.949249554E+01’ 
351, 0.101592350E+02, 
352, 0.108298464E+02/ 
353, 0.5949458313E+01, 
354, 0.101612844E+02, 
355, 0.108299999E+02; 
356, 0:949667072E+01, 
357, 0.101633348E+02, 
358, 0.108299999F+02/ 
359, 0.114600000E+02/ 
360, 0.115368347E+02, 
361, 0.442086369E+00, 
362, 0.442086369E+00, 
363, 0.442086369E+00; 
364, 0.442086369E+00, 
365, 0.442086369E+00, 
3, SABES. 

*NSET, NSET=MERIDLL, GENERATE 

3,,255,9 5. 

ASET, NSET=MERID22 , GENERATE 

4,256,9 E 

ASET; NSET=MERID34 , GENERATE 

PSE’ x S cneae 
/NSET=MERID4 

6,258 ,,9 = 

ASET, NSET=MERIDS6 , GENERATE 

leeeere a CHUTE 
NSET=MERID6 

8 260,9 5 

NSET, NSET=MERID79 , GENERATE 

9,261,9 : 

*NSET/ NSET=MERID 

ANSET 

2120,255,9 

*NSET, NSET=STRIP22 , GENERATE 

Saupe ee-ctoreya eens 

212,257,9 ’ 
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sNSET NSET-STRIPSS, ,GENERATE 
eS unite beers 
sNSETN 9 oe 
*NSET, NSET=STRIP79 , GENERATE 
218,261 9 
smi STRIP22,STRIP34 , STRIP45, STRIPS6, STRIP67 , STRIP79 


308 cen 313,355, 358, 360 
ey NSEISSHIZO, GENERATE 


NSET, NSETSCORE 

265,266,271 272 

* 

1,Z0,MERID11 , MERID22 ,MERID34 ,MERID45 
TYPES=C 


0. 2a ith 542023210 10. : 
TORPELD, 1 
323) 3 


*ELEMENT, TYPE=STRI35 _, ELSET=SHLP15 


A 3, 2 1 
2; 4; 3; I 
ci 5) 4! ir 
4) 6, 5) 1 
5; 7; 6, I 
6, 8) 7, I 
7, 9, 8, 1 
8 10 9 lL 
* , TYPE=C3D4 > =BIADSLD 
303, 1, 265, 263, 361 
304; 1, 265, 361, 362 
305, 1, 265, 362, 363 
306, 1, 265, 363, 364 
307; 1, 265, 364, 365 
308/ 1, 265, 365, 365 
309/ 1, 265, 366, 367 
310; 1, 265, 367, 314 
311, 1, 361, 263, 3 
312; 1) 263, 2; 3 
313, 1) 362, 361, 4 
3147 1) 361) 3) 4 
315; 1) 363) 362) 5 
F-9 


316, 1, 362, 4, 
317, {i 3647 363; 
je meet 
320, 1) 365) 6, 
321, 1) 366, 365, 
33,388) seh 
324; i? 3367, 8) 
ee be Sey 3ege 

* TYPESSTRI35 —’ 
233, 263, 1 2 
235, 265, 1) 263 
a ae 

*ELEMENT, TYPE=S4R5 : 

Gg, 2, 3; 2 a2, 
10) a 4. 13) 
ey FY 2! 1 
, , 7 a 
13) 6, 7, 16) 
14) 7 oS a7) 
te’ 8'_—49! 18! 
*ELEMENT, TYPE=S485 : 
5 y sealed Saari be Pea 
13) .. hae. og) 
7 Fd a tA 
21, 157 6, 25) 
22, 16, i7; 26, 
23,0 -T7,> Be ~27) 
24, 18, 19, 28, 
35; 20, 21, 30, 
D6... “21, 22, Bl, 
O7... 22; 223, . 32, 
28, 23, 24) 33; 
39’ = 34 35) 34), 
30, 25, 26, 35, 
a a a 
* TYPE-S4R5 |, 
, , 
33, 29, 30; -39, 
34’ 30 31) 40) 
35°. 31, 32-41, 
36°. 82, a3," 42, 
37, 33; 34) 43, 
38/ 34) 357 44) 
39' 35, 36, 45; 
40; 36, 37; 46, 
41, 38, 39, 48, 
42; 39; 40, 49, 
43, 40, <All > 50, 
44  4y) 43% BT 
45, 42, 43) 52, 
46, 43; 44, 53, 
47, 44, 45, 54) 
48° 45, 46, 55, 
49° 47) 48/ 57; 
50, 4g 49° 58) 
51; 49; 50, 59; 
52, 50, 51; 60, 
53/51; 52,61 
54, 52, 53, 62, 
55; 53) 4654) 63, 
F-10 
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13, 
308! 14 
Tso’ 304; 203, 215 
; 205; 6 
Tou! 206, 207, 218, 
to3! 207, 208; aids 
193! 209; 210, ge, 
Toa 210, 211; 30, 
Tos! 211, 212) 3a, 
T3e! 212; 213, 3381 
Tos’ 313 214, 22 
198! 2145-215, 334, 
Tos! 215, 216, 333) 
500’ te 217, 336, 
309" 218; 219) sar 
303° 219; 220) 3231 
303! 220, 221; 330, 
304! 221, 222, 33h, 
308! 222, 223, $34" 
30e" 223, 224, gai) 
309! 224, 225; 334, 
308" 225, 226, gan" 
309" 227, 228, rere 
398" 228, 229) gag! 
ant 229, 230) a38) 
212 1 33h, 241 
34: 2 334" 242, 
t 

215, 23, 333! 244) 
; 234) 8 
315" 236, 237, 246) 
318" 237, 238, ga, 
318" 238, 239, 348) 
330! 239, 240) 343) 
330" S40! 241, 3201 
333! 241, 242) 3ph, 
353" 242, 243; guar 
334 343! 244) guar 
3340 245, 246, 332) 
33e' 34! 247, 26 
ssc $45" 248, a4) 
338! 248; 249) 338, 
338" 249; 250) 328" 
330" 250, 2517 20, 

? ; 252, 
231; 323! 253 262, 
; 336 263! 264; 206, 
, 263/ a 
236, 263, 26 ; 388" 
338" 267, 268, 2701 
339" 264, 269, aint 
340’ 5e9! 270, 321 
349! 268, 273) 387 
343! 273, 274; 34, 
343! 274, 275, guar 
344" 275, 276, 2g6, 
348' 276, 277, ee 
342! 399° 8, 88, 
349! 278, 279, coy 
348! 279, 280, gar 
349° 280, 281, agi, 
350! 281, 282, 333, 
3a" 270, 283, 331 
333! 283, 284; a34, 
$23" 284; 285, 238) 

A ; 286 

323" 285, : 
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MINNINA 
SERCERRS 


221 


255, 286, 287, 297, 
256, 287, 288, 298) 
257, 288, 289, 299, 
258, 289, 290, 300, 
259, 290, 291, 301, 
260, 291, 292, 302, 
261, 282, 303, 306, 
262, 303, 304, 307, 
263, 304, 305, 308, 
264, 292, 306, 309, 
265, 306, 307, 310, 
266, 307, 308, 311, 
267, 305, 312, 313, 
269, 10, (19, 315, 
271, 314, 315, 266, 
272, _19, 316, 318, 
273, 316, 317, 319) 
274; 315, 318, 272, 
275, 318, 319; 272, 
276, 317, 320, 330, 
277, 320, 321, 331, 
278, 321, 322, 332, 
279, 322, 323, 333, 
280, 323, 324, 334; 
381; 324, 325, 335, 
282, 325, 326, 336, 
283, 326, 327, 337, 
284, 327, 328, 338, 
385, 328, 329, 339) 
286, 319, 330, 340; 
287, 330, 331, 341; 
288, 331, 332, 342) 
289, 332, 333, 343; 
290, 333, 334, 344) 
291, 334, 335, 345, 
292, 335, 336, 346; 
293, 336, 337, 347, 
294; 337, 338, 348, 
295, 338, 339; 349; 
296, 329, 350, 353; 
297, 350, 351, 354; 
298, 351, 352, 355) 
399/ 339, 353, 356, 
300; 353, 354, 357, 
302’ 380’ 388 360! 

*ELEMENT, TYPE=SPRING2 | 
327, 282, 119 
328, 303, 128 
329; 304, 13 
Be ig. tee 

*ELEMENT 
332, 329, 127 
333; 350, 136 
334; 351, 145 
33’ 386’ 1e3 

6 

*ELEMENT, TYPE=SPRING2 , 
337, 282, 119 
338; 303, 128 
339/ 304) 137 
340; 305; 146 
341), 312, 155 
342/ 329, 127 
343; 350, 136 
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SSPRING, ,ELSET=YSPRG 

1. F06 

*SPRING , ELSET=2SPRG 

1! E06 

¥SPRING, ELSET=XSPRG 
1/02 

*MATERTAL, NAME=FIBERTTE 
2.4E+06,0.28,72. 

*DENS 


1.682F=-4 
. , NAME=H250 


Cc 
5.51E+04,0.1,72. 
*DENSITY 
2. 345-05 


1S 


TRAMEOETOS, TIME=D, VALUE=A 
-450E-05" 9. OOOE+00’ 1.770E~04-3 . 728E+02 
- 197E-03-1.033E+02 aecoeee 3-1. 006E+02 
715E03-5, 6E+01 3 aaa eee 402E+01 
26E- 93-9. 194F+01 4:734E-03-9:121E+0 


33 

8 . E+) 
ee ~03- -8:580E+01 8:212F-03-8.4 aaoroe 
«1 


oO 
Sa 
5 
te 
coos 
to- 
K@ 
be 
le) 
iu 
Dw 
iE 
dob 
200} 
Bee 
— 


D3E-03-8.223E+01 1.026E-02-8.158E+0 
77E-02-8.031E+0L 
3, TIME=D, VALUE=A 
450E-05'0.000E+00'1.770E-04 0.000E+00 
S7E~03-7 899E+01 1:562E 1 a ee BE 
375E+0 ~03-7.304E+01 


i 
1. 
5E+01 3.113E 
3-7-1970) 4. 7348-03-74 6E+O1 
8. 
i: 


ba 


Z0gR 0376 BRE tOt 
212E-03-6.608E+01 
026E-02~6. SBeE TOL 


M00 
NES 
ree 
d 
WU) 
og had 
WwW 
» 
fa 
& 
| 


TE~02 


,NAME=EIG4, 
-450E-05"0.000E+00’ 1. 
9TE-03-7.900E+01 1. 
15E-03 LE+0 i 


x 
e 
od 
peu 
bel 
db5¢ 
Pr 
Peo. 
22 
dh 
Sto 
Bs 
ze 
i 


-5.78 12 
» {538 -03"3- 445E+01 1. 626E-02-8. AB3ETOL 
TE 5.261 E+01 


 NAME=ELGS, TIME=D, VALUE=A 

50E-05’0.000E+00'1.770E~04 0.000E+00 
7E-O3-1.110E+O2 ].562E-03-6.642E+01 

3.113E-03-5.364E+0) 
26E-03-5.295E+01 4.734E-03-5.256E+01 
DE-03- 6.264E-03-5. 018E+01 
8.212E-03-4.744E+01 
1 1.026E-02-4.413E+01 


HIN JJ~100 Lo apt 
0 ~] N 


5 
6. 
2 ee ee 1.141 


5.040E-04-1.637E+02 


NNN 


04 . 
938E-03-9 
515E 


5.040E-04 0.000E+00 
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. 2 
asap esos seaeet 


3 
BE-03-4  559E+01 
127E-02-4.324E+01 


8.440E-04 0.000E+00 


a meee Le LULU 


1.197E-03 0.Q00E+00 1.562E-03-9.336E+01 1.935E-03-5.553E+01 
3° 715E-03-4.615E+O1 3. 113E-03-4.535E+01 3.515E-03-4.497E+01 
4.326E-03-4. qB7ErO1 4: 734E-03-4.452E+01 5:116F-03-4.377E+01 
5. 895F-03-4.256E+01 6.264E-03-4.2 1 6.741E-03-4 .157E+0 
7° 713E-03-4.024E+01 8:212E-03-3.944E+01 8:722E-03-2 .853E+01 
92753 3E-03-3. 669E+01 1.026E-02-3.604E+O1 1:°76E-02-3 .554E+01 
1 177E-02-3.481F+01 
x TORE, Nave 

4.450 0.000E+00' 1. SF OBLOd 0.000E+00 5.040E-04 0.000E+00 
1,197E 3 8° 00100 17562E-03 0.000E+00 1:938E-03-7.861E+01 
4-394 93-4 0a les L 32113E-03-3.889E+01 3:515E-03-3-846E+01 
4..326E-03-3.849E+01 4:734E-03-3.817E+01 5. 116E-02-3.737E+0L 
5.895E-03-3.611E+01 6.264E-03-3.572E+02 6. JAE 03-3. S14EtOL 
7. 713E-03-3.380E+01 8.212E-03-3.298E+01 8:722E-03-3:203E+0 
9.753E-03-3.013E+01 1.026E-0z-2.946E+01 9: O52B-02 3 AOaEtOL 
gdh] TE-02-2, 820E+01 

rere iise 0.000E+00' 1.770E-04 0.000E+00 5.040F-04 0.000E+90 
L197E-03 0:000E+00 1.562F-03 0.000E+00 1:938F-03 0: 000EHI0 
3° 715E-03-3 97 LEFO1L 3. 113E-03-3.437E+01 3.515F-03-3. 355E+01 
4.326E-03-3.364E+01 4:734E-03-3.3 1EHOL 5: ~03-3.24°E+01 
51895E-03-3. 107E+01] 6.264E-03-3.06 6. 741E-03-3 . OO7E+O1L 
**713E-03-2-.869F+01 8:212F-03-2. OS4E TOL 8. 722F-03-2 . 685E+01 
9: 733F-03-2 -485E+01 1 026E-02-2.415E+01 1:076E-02-2.362E+0L 
shiek? TE-O2-2 , 285E+0 4 

4.450E-05'0. taster! 1.770E-04 0.000E+00 5.040E-04 0.000E--00 
1.197E-03 0.000E+00 1.562E-03 0.000E+00 1:938E-03 0.000E+00 
2° FISE-03-5.734E+0L 3. 113E-03-3.452E+O1 3.515E-03-3.023F+01 
4.326E-03-3.007E+01 4:734F-03-2.969E+01 5°116F-03-2:865E+01 
5 BO5E-03-2.715E+01 6.264E-03-2.672E+01 6.741E-03-2 . 6CSE+01 
7°713E-03-2.462E+01 8.212F-03-2.372 8.722E-03-2 . 269F+01 
9:753E-03-2.058E+01 1.026E-02-1.984E+01 1:076E-02-1.929E+01 
1.177E-02-1.849E+01 
* E, 10 , TIME=D, VALUE=A 

4.450F-05'0.000F+00 1.770F-04 0.000E+00 5.040F-04 0.000E+00 
1:197E-03 0:000F+00 1:562E-03 0.000E+00 1:938E-03 0:000E+00 
3° 715E-03 0.000E+00 3.113E-03-5.031E+01 3.515E-03-3.094E+0 

4. 326E-03-2.760E+01 4:734E-03-2.710E+01 5:116E-03-2:583E+01 
5.895E-03-2.410E+01 6.264E-03-2.362E+01 5.741E-03-2.293E+01 
7°713E-03-2.135E+01 8.212E-03-2.039E+01 8.722F-03-1.529F+01 
92753E-03-1.705E+O1 1.026E-02-1.627E+0l 1.076F-02-1.568E+01 
adie! ZEO2-1. 485 tol, 

4.450E-05'0.000E+CO 1.770E-04 0.000E+00 5.040E-04 0.000E+00 
1.197E-03 0:000E+00 1:562E-03 0.000F+00 1:938E-03 0:000E+00 
3°715E-03 0.000F+00 3.113E-03 0:000E+00 3.515F-03-4 /526E+01 
4.326E-03-2.614E+01 41734E-03-2.530E+01 5.116E-03-2:365E+0 
5.895E-03-2. 155E+01 6.2645 -03~2. 1O4ETOL 6.7422~-03-2.025EH01 
FIFIZE-A? Tos4E+O1 8:212F-03-1.750E+01 8:7226-03-1.632E+01 
9-753E-03-1. 393E+0 : 1:026E--02-1.309E+01 1.076F-02-1.248E+01 
x LE, TIME=D, VALITESA 

4.4506-05 0.000E+00 1.770E-04 0.OUCE+00 5.040E-04 0. QQ0E+00 
1:19 3 O000E+00 1:562E-03 0:000E+00 1:938E-03 0:000E+0 
2° 715E-03 0:000E+00 3.113E-63 0.000E+00 3.515E-03 0: OUsEO0 
4.326E-~03-2. +01 4:734E-03-2.404E+01 5: 116F-03-2: 165E+0. 
5° 895E-03-1.908E+01 6.264E-03-1.851F+01 6.741E-03-1.767E+O] 
72713E-03-1.576" +01 8.212F-03-1.462E+01 8:722F-03-1.333E+01 
9: {53E-O3-1- O77L+01 1:026E-02~9.883E+00 1:076E-02-9. 230E+00 
13177E-02-8 : 316E+00 
x SDE, NAME=ETG13 =D, VAIUE=A 

4.450E-05" Q. OO0E+00 1.776E-04 0.000E+00 5.040E-04 9. Q00E+00 
1,197E-03 0:090E+00 1.562E-03 0.000E+00 1:938E-03 0:000E+00 
27 715E-03 0.CCOF+00 3:113E-03 0.0C0E+00 3.515F-03 0. 000EHO0 
4.326E-03-3.981E+01 4.734E-03~2.474E+O1 5.1I16E-03-1.907E+01L 
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2 
48E~03-3.101E+01 


Cis ade 


2 
127E~-02-1.885E+0" 


“1. 808E+0). 


455; 
38b 0374 . SoH 


TE-02~8 72> 5+00 


E-04 0.000E+00 
E-03 0.QQ0E+UC 
me 0. Q00F+00 

3-1.694E+01 


5.895E-03-1.586E+01 6.264E-03-1.526F+01 §.741E-03-1. 436E+01 7.223E-03-1.336E+01 
7°713E-03-1.228F+01 8:212F-03-1. 104E+01 8.722F-03-9.651EH00 9.248E-03-8. 148E+00 
9.753E-03-6.881F+00 1.026F-02-5.938E+00 1.076E-02-5.248E+00 1:127E-02-4.714E+00 
1: 177E-02-4.290E+00 
* FIG14,T 
4.450F-05' 0.000E+00 1.770E-04 0. QOQEHO0 5.040E-04 0.000E+00 8.440E-04 0.000E+00 
1:19 0:000EF+00 1:562E-03 0:000E+00 1.938E-03 0:000F+00 2:323E-03 0:000E+00 
3.7 0: Q000E+00 3:113E-03 Mestostiets 3. 03 0.000E+00 3.920E-03 0:000E+00 
4.326E-03 0.000E+00 4:734E-03-3.402F4+01 5.116F-03-1-449E+01 5:504E-03-1.040E+01 
5.89 ~9.370F+00 6.264E-03-8.964E+00 8. JE-03-8.144E+00 7.223E-03-7.108E+00 
7°713E-03-5.989E+00 8.212F-03-4.679E+00 8.722F-03-3.191F+00 9.248E-03-1.584E+00 
9.753E-03-2.249F-0l 1:026E-02 0: 000E+00 1.076E-02 0:000E+00 1:127E-02 0.000E+00 
Rey 502 0.0 E+O6 
4.4506-05 0.000E+00 1.770E-04 0.000E+00 5.040E-04 0.000EH00 8.440E-04 0.000F+00 
1.197E-03 0.000F+00 1.562E-03 0:000F+00 1.938E-03 0.000E+00 2.323E-03 0.000E+00 
2:715E-03 0.000F+00 3:113E 0. 000E+00 3. ~03 0:000E+00 3.920E-03 0.000E+00 
4.326E-03 0.000E+00 4:734E 0.000E+00 5.116E-03-1.658E+01 5:504E-03-4.784E+00 
5.895E-03-2.315F+00 6.264E-03-2.282F+00 6.741F-03-1.703E+00 7.223E-03-6.985E-O1 
7.713E-03 0.000E+00 8:212F-03 0.000E+00 8.722F-03 0.000EH00 9:248E-03 0.000E+00 
9.753E~03 0.000E+00 1.026E-O02 0.000E+00 1.076E-02 0.000E+00 1.127E-02 0.000E+00 
1:177E-02 0.000E+00 
* 16 VAIUE=A 
4.4506-05'0.000E+00 1.770E-14 0.000E+00 5.040E-04 0.000E+00 8.440E~04 0.000F+00 
1: 03 0.000E+00 1.562E-03 0. 000F+00 1.938F-03 0:000F+00 2:323E-03 0.000E+00 
3° 715E-03 0:000E+00 3:113E-03 0: 000E+00 3.515F-03 0-.000EH00 3.920E-03 0.000E+00 
4.326E-03 0.000E+00 4.734E-03 0.000E+00 5.116F-03 0:000E+00 5.504E-03-1.493E+01 
5 .895E-03-4.452F+00 6.264E-03-3.184E+00 6.741F-03-2.419F+00 7.223E-03-1. 102E+00 
7°713E-03 0.000E+00 8:212F-03 0:000EH00 8:722F-03 0.000E+00 9:248F-02 0.000E+00 
9.753E-03 0.000E+00 1:026E-02 0.000E+00 1:076E-02 0:000E+00 1:127E~02 0:000E+00 
1:177E-02 0.000E+00 
* ,NAME=ETG17 
4.450E-05'0.000E+00 1.770E-64 0.000E+00 5.040E-04 0.000E+00 8.440E-04 0.000E+00 
1.197E-03 0.000F+00 1.56 3 0.000E+00 1.938E-03 0: 000E+00 2:323E-03 0.000E+00 
2°715E-03 0.000F+00 3:113E-03 0:000E+00 3.515F-03 0:000F+00 3.920F-03 0.000E+00 
4.326E-03 0.000F+00 4:734E~03 0: 000E+00 5.1126F-03 0:000E+00 5:504E-03 0.000E+00 
5.895E-03-1.207E+01 6.264E-03-6.004E+00 6.741F-03~3.419F400 7.223F-03-1.489E+00 
7.713E-03 0.000E+00 8:212F-03 0.000E+00 8.722E-03 0:000E+00 9.248E-03 0.000E+00 
9:753E~03 0.000F+00 1:026E-02 0: 000EH)0 1.076E-02 0:000E+00 1:127E-02 0: 000E+00 
1.177E-02 0.000F+00 
* ; 18 , TIME=D 
4.450F-05'0.000F+00 1.770E-64 0.000E+00 5.040E-04 0.000E+00 8.440E-04 0.000E+00 
1:197E-03 0:000E+00 1.562E-03 0:000E+00 1:938F-03 0:000E+00 2.323E-03 0.000E+00 
3° 715E-03 0.090E+00 3-.113E-03 0. 000E+00 3.515F-03 0:000F+00 3.920E-03 0:000E+00 
4.326E-03 0.000E+00 4.734E-03 0.000E+00 5.116F-03 0:000F+00 5:504E-03 0.000E+00 
5.895E-03 0.000E+00 6.264F-03-1.798E+01 6.741F-03-5.583E+0 3-22 23E-03-9.063E-01 
7.713E-03 0.000E+00 8:212F-03 0.000E+00 8.722F-03 0:000F+00 9:248E-03 0. 000E+00 
9.753E-03 0.000F+00 1:026E-02 0: 000EH00 1:076E-02 0:000E+00 1:13 7E-02 0. 000E+00 
ard? E702 9. ,000E+00, 
4.450E-05 0.000E+00 1.770E-64 0.000E+00 5.040E-04 0.000F+00 8.440E-04 0.000E+00 
1.197E-03 0.000F+00 1.562E-03 0.000F+00 1.938E-03 0. 000E+00 2.323E-03 0.000E+00 
2:715E-03 0.000E+00 3.113E-03 0:000E+00 3.515F-03 0.000E+00 3.920E-03 0.000E+00 
4.326F-03 0.000F+00 4:734E-03 0: 000F+00 5.116F-03 0:000F+00 5:504E-03 0:000E+00 
5°895E-03 0.000E+00 6.264F-03 0.000E+00 6.741F-03-1.021E4+01 7.223E-03 0.000E+00 
7.713E-03 0.000E+00 8.212F-03 0.000E+00 8.722F-03 0:000E+00 9.248E-03 0. 000E+00 
9.753E-03 0.000E+00 1.026E~02 0: 000E+00 1. ~02 0.000EH00 1.127E-02 0: 000E+00 
1.177E-02 0.000E+00 
* ,NAME=ETG20 , TIME=D, 
4.450E-05'0.000F+00 1.770E-04 0.000E+00 5.040E-04 0.000F+00 8.440E-04 0.000E+00 
1.197E-03 0.000F+00 1.562E-03 0:000E+00 1:938E-03 0:000E+00 2:323E-03 0.000E+00 
2:715E-03 0.000E+00 3:113E-03 0:000E+00 3.5 0.000E+00 3:920E-03 0:000E+00 
4.326E-03 0.000F+00 4:734E-03 0:000F+00 5.116E-03 0:000E+00 5:504E-03 0.000E+00 
5.895E-03 0.000E+00 6.264E-03 0.000E+00 6.741E-03 0:000E+00 7.223E-03-8.116E+00 
7°713E-03 O:000F+00 8:212E-03 0:000EH00 8:722F-03 0:000E+00 9:248E-03 0. 000E+00 
9.753E-03 0.Q00E+00 1:026F-02 0:000EH00 1:076E-02 0:000E+00 1:127E-02 0:000E100 
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a dZ7En02 0.000E+00 hed ; 
STEP, AMPLITUDE=RAMP, LINEAR, INC= CLE= 
TRANSTENT RESPONSE OF’NOSECAP’TO IST WATER-ENTRY STEP PRESSURE 


*DYNAMIC 
4.45E-05 4.45E-05 


DLOAD, AMPLITUDE=ELG1 
EIG1 PA. 


*STEP, AMPLITUDE=RAMP, LINFAR, INC= 1, CYCIFE 1 
TRANSTENT RESPONSE OF’ NOSECAP’TO 2ND WATER-ENTRY STEP PRESSURE 


*DYNAMIC 
1.325E-04,1.325E-04 
bal FIG1 


DLOAD, AMPLITUDE= 
EIGI,P, 1. 
FLED Pil 
*END STEP 


*STEP, AMPLITUDE=RAMP, LINEAR, = 1, CYC #1 

TRANS TENT RESPONSE OF’ NOSECAP TO 3RD WATER-ENTRY STEP PRESSURE 
3.270E-04 ,3.270E-04 

*DIOAD, AMPLITUDE=EIG1 


D 
EIGI,P/1. 
*DLOAD) AMPLITUDESELG2 


*STEP, AMPLITUDE=RAMP, LINEAR, INC= 1, CYCIF= 1 
TRANSTENT RESPONSE OF’ NOSECAP’ TO 4TH WATER-ENTRY STEP PRESSURE 


DYNAMI 
3.400E-04 ,3.4000E~04 
*DLOAD ELG1 


*STEP, AMPLITUDE=RAMP, LINEAR, INC= 1, CYCIE= 1 
TRANS TENT RESPONSE OF’ NOSECAP’TO 5TH WATER-ENIRY STEP PRESSURE 


ere) 04 3.533E-04 
*DLOAD, AMPLITUDE=ELGL 
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*END STEP 
*STEP 


a eyhseehh 3.651E-04 
*DLOAD, AMPLITUDE=FIGL 


ee : eEnO4 3.756E-04 

*DLOAD, AMPLITUDE=FIGL 
ELG1,P,1. 
*DLOAD) AMPLITUDE=ELG2 


oie bids 


*STEP, AMPLITUDE=RAMP, LINFAR, INC= 
TENT RESPONSE OF’ NOSECAP’ TO 8TH 


TRANS 
SS. 3.846E-04 
*DLOAD, AMPLITUDE=ELG1 
FIG1,P,1. 


*STEP , LINEAR, INC= 
TRANSTENT RESPONSE OF’ NOSECAP’ TO STH 


3 EOpioes 4,3.922E-04 
*DLOAD, AMPLITUDE=EIGL 
FIG1,P,1. 

*DLOAD, AMPLITUDE=ELG2 
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AMPLITUDE=RAMP, LINEAR, INC= 
TRANSLENT RESPONSE OF’ NOSECAP’TO 61TH 


1, cy 


WATER ENTRY it grip PRESSURE 


1, CVF 1 
WATER-ENTRY STEP PRESSURE 


1, CYCIE= 1 
WATER-ENIRY STEP PRESSURE 


END’ STEP 

*STEP, AMPLTTUDESRAMP, LINFAR, INC= 1, CYCIF= 1 

TRANS TENT RESPONSE OF’ NOSECAP’TO IOTH WATER-ENIRY STEP PRESSURE 
3, 980-04 3.980E-04 


DLOAD, AMPLITUDE=ELGI 
FIG1, P71 


* 

4.023E-04,4.023E-04 
*DLOAD, AMPLITUDE=ELG1 
FIGL,P,1l. 
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*END STEP 
*STEP, AMPLITUDE=RAMP, LINEAR = 1, CCE 1 
TRANSIENT RESPONSE OF’ NOSHCAP’TO 12TH WATER-ENTRY STEP PRESSURE 
4. Oh Oe O4 4.049E-04 


DLOAD, AMPLITUDE=EIG1 
FIGL,P/1. 


*STEP, AMPLITUDE=RAMP 


4 O57E-04, 4 057E-04 

4037 : é 
Bere ne 
*DLOAD) AMPLITUDE=EIG2 


oto as 4,4.078E-04 
*DLOAD, AMPLITUDE=FIG1 


INC= 


LINEAR, INC= 
TRANS TENT RESPONSE OF’NOSECAP’TO 13TH WATER-ENIRY STEP PRESSURE 


*STEP, AMPLITUDE=RAMP, LINEAR, INC= 
TRANSTENT RESPONSE OF’ NOSECAP’TO 14TH 
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STEP PRESSURE 


HGLP 1 

Bret Br 1s MPLITUDEX FIG3 

ELG3,P,1. 

*DLOAD) AMPLTTUDE=ELG4 

Thr aab! 1s MPLTTUDE=ELGS 

FIGS, P,1. 

*DLOAD, AMPLITUDE=ELG6 

GEeCA 1s MPLTTUDE=EIG7 
s 

TOT AAD’ SAMPLTTUDE=ELGS 
t 

EOPeADY 1s MPLTTUDE-ELG9 

BECO (Pr Ls sor TUDEXETGIO 

Eine PL prTTUDE= FIGIL 

FE Pod or ESIC? 

Fee Pr or TTUDESEIGL 

EIG13,P,1. 

*DLOAD, ’ AMPLITUDESEIG14 

EIGI4, PL. 


*STEP, AMPLITUDE=RAMP, LINEAR, INC= 2, CYCIF= 1 
TRANSTENT RESPONSE OF’ NOSECAP’TO 15TH WATER-ENTRY STEP PRESSURE 


SOR ORLo4 3.829F-04 
*DLOAD, AMPLITUDE=EIGL 
FIG1,P,1l. 

*DLOAD) —AMPLITUDE=EIG2 


*STEP, AMPLITUDR=RAMP, LINEAR, INC= 1, CYCFE= 1 
"NOSECAP’TO I6TH WATER-ENTRY STEP PRESSURE 


F-23 


3.877E-04,3.877E-04 
*DIOAD, DE=ETG 
EIG1,P,1. 
*DLOAD, AMPLITUDE=ELG2 
FIG2,P)1. 
*DLOAD, AMPLITUDE=EIG3 


* 
*STEP, AMPLITUDE=RAMP, LINEAR, INC= 1, CYCIER= 
TRANSTENT RESPONSE OF NOSECAP’TO 17TH WATEXCENTRY step PRESSURE 


* 
3.911E-04,3.911E-04 
*DLOAD, AMPLITUDE=ELG1 


*STEP, AMPLITUDE=RAMP, LINFAR, INC= 1, CYCIF= 1 
TRANS TENT RESPONSE OF NOSECAP’TO 18TH WATER-ENTRY STEP PRESSURE 


DYNAMIC 
3.690E-04, 3 .690E-04 
*DLOAD FIGL 
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USER INSTRUCTIONS— WATER 
ENTRY STRUCTURAL TECHNIQUE 


APPENDIX G 


USER INSTRUCTIONS-WATER ENTRY STRUCTURAL TECHNIQUE 


These user instructions for WEST are written assuming that the 
user is familiar with both PATRAN and ABAQUS. Therefore, no effort is 
expended on their use in these instructions. The user instructions for 
WEST are divided into three parts: 

Part 1 provides the instructions for the generic use of the 
ENTRY code, without any pre~ or post-processing nor finite element 
analysis of the structure. These instructions are extracted from the 
original users manual for ENTRY, which is contained in reference 19. 

Part 2 lists the special instructions a user will need to generate 
a finite element mesh that can accommodate or bypass the restrictions 
imposed by the ENTRY code. 

Part 3 provides instructions on how to run the translators PATENTR 
and ENTPRES, and how to activate the added output options for element 
pressure- and load-time history files in an ENTRY run. 


PART 1 


INSTRUCTIONS FOR THE GENERIC USE OF THE ENTRY CODE 


The current version of the ENTRY code can be applied to arbitrary 
bodies. The grid describing the entry body may contain up to 750 nodes 
and 500 elements. However, execution will terminate when more than 300 
of these elements become submerged. These instructions describe the 
available program options, necessary input cards and output format. 


Program Options and Required Input 


Program input can be divided into three parts. In the first, the 
basic program options are specified: 


Card No. Variable Format 
1 CONSTANT or VARIABLE body orientation 2A4 
2 PRINT or DON'T PRINT 3A4 
3 ASYMMETRIC or SYMMETRIC mode 3A4 


Under the CONSTANT body orientation option the entry model is 
assumed to retain its initial orientation and velocity throughout the 
entry process. The natural problem variable in this case is depth 
rather than time. With little increase in computational time, 
pressures and forces can be evaluated for a number of different wetting 
factors, Cy. The VARIABLE body orientation option allows the velocity, 
orientation, wetting factor and time increments between steps to be 
varied. The only restriction is that the angular velocity of the body 
must be small enough to insure that the depth of the body increases 
monotonically in time. The maximum number of stevs is limited to 49. 


The PRINT option is used to obtain flow field and element 
information at each step of the calculation. It is applied only for 
diagnostic purposes. The second option, DON'T PRINT, is recommended 
and produces only grid information and the final pressures and forces 
on the model. 


If the SYMMETRIC mode option is used, the entry model is assumed 
to possess planar symmetry about the y-z plane. The ASYMMETRIC option 
does not assume any symmetry and hence can be applied to arbitrary 
bodies. This mode is also used on symmetric bodies where Vy is non- 
zero. 


The second set of input cards describes the entry conditions. The 
required information differs depending on whether the CONSTANT or 
VARIABLE body orientation option is used. For the CONSTANT option the 
following data cards are required: 


Card No. Variable Format 
4 IMAX, D, VENTRY, ANG, SUMT, HMIN, DH, ALPHA I5,5X,7F10.0 
5 CGL, FCF, ANGB, NNLD 3F10.0, I5 
6 NCW, CW(1), CW(2).... (CW(NCW) ) I5,5X,7F10.0/ 
(8F10.0) 
7 omit 


These variables are defined as follows: 


IMAX Number of steps at which pressures and loads are calculated. 
The present calculative procedure inserts the model into the 
water in a series of steps, each at a greater depth than the 
preceding one. When the step count becomes greater than 
IMAX, execution is terminated. 


D Diameter (in inches). This quantity is only used for 
calculating force coefficients. 

VENTRY Entry velocity in (in./sec.). 

ANG Orientation of the model (in degrees) relative to the water 


surface (see Figure G-l). 


SUMT Program time limit. This variable is not used in the VAX 
implementation of this code. Set equal to zero. 


HMIN Initial body depth (i.e., measured from the lowest point on 
the body). This parameter is zero if the loads are 
calculated from the time of initial wetting. Note that if 
this variable is not zero pressures and forces are first 
calculated at HMIN+2DH. This parameter allows pressures and 
forces at a particular depth to be determined without 
calculating the entire force-time history from initial 
wetting. 


DH Increment in depth (in inches) between successive steps. It 
is necessary to coordinate this variable with the specified 
model grid which is defined on the last set of data cards. 
The following apply to determining DH: 


a. OBLIQUE ENTRY WITH STANDARD GRID OPTION. DH should be 
picked so that the average element is submerged in two steps. 
On models of complex shape this criteria can only be 
satisfied in the mean and primary consideration should be 
given to the portion of the body which experiences the 
greatest load. Generally this will be on elements whose 
plane is perpendicular to the direction of motion. 


b. VERTICAL ENTRY WITH STANDARD GRID OPTION OR OBLIQUE ENTRY 
WITH THE OGIVE OPTION. For vertical entry or if the OGIVE 
grid option is used, elements will have a pair of sides 
parallel to the water surface. In this case it is important 
to choose the step size very precisely so that each element 
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FIG.G-1 TERMS DEFINED 


EFFECTIVE 
PLANNAR 
SURFACE 


m —> 


FIG. G-2 PROFILE OF CONE GRID 


ANG 


ALPHA 


CGL 


FCF 


ANGB 


NNLD 


NCW 


CW 


will be submerged in exactly two steps. To insure that the 
top row of elements is included in an unmodified state in the 
code and that the next row of elements is excluded, the 
actual water surface should fall a small distance « above 
the upper edge of the top row of elements to be included as 
shown in Figure G-2. Here 0 <« < Ah where Ah is defined 


by 
th-= [average element area 
1000 


Angle of attack (in degrees). See Figure G-l. 
z' coordinate of the center of gravity. See Figure G-1l. 


Pressure correction factor on elements with a modification 
code of 1. For the oblique entry of blunt bodies (nose 
length/diameter < 1) set to unity. For other cases use a 
value of 0.67. 


Yaw angle (in degrees). Velocity components in the x, y, z 
directions are Vysin (ANGB), -Vycos(ANG)cos(ANG+ALPHA) and 
-Vzcos (ANG) sin(ANG+ALPHA) . 


Number of no-load elements. 


Number of wetting factors to be used. Since the most 
appropriate value may not be clear, for little extra 
computational cost, pressure and loads may be calculated for 
several different wetting factor values. 


Wetting Factor. This parameter describes the rate of surface 
rise and is equal to the ratio of h/h' defined in Figure G-3. 
For best results, the test cases reported on in Reference 19 
should be used as a guide. An approximate rule for 
determining this parameter is as follows: 
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(1) POINTED BODIES (ALSO INCLUDES SLIGHTLY BLUNTED ONES). 
Determine the angle, %. , between the tangent to the body 
surface and the body axis at both the nosetip and base of the 
nose. At the nose neglect any effect due to body blunting. 
Insert the two resulting values of in radians into: 


(1 -.3966. +.1870, 2 -.1240¢ 3) 


Average the two calculated values of Cy to obtain the final 
one to be used in the code. If ALPHA is non-zero increment 6, 
by ALPHA. 


(2) FLAT PLATES. Use a value of 1.45 for ANG>45 degrees and 
1.55 for ANG<45 degrees. 


(3) SPHERICAL BODIES. Use a value of 1.55 for near vertical 
entry and 1.35 for oblique entry. 


The classification of an arbitrary body into one of the above 
categories is a matter of experience. On complex shapes 
classification should be based on the portion of the body 
sustaining the majority of the impact loading. 


If the VARIABLE body orientation option is used the following data 
cards are required: 


Card No. 


NOS 


Variables Format 
IMAX, D, VENTRY, ANG, SUMT, HMIN I5, 5X, 7F10.0 
CG., FCF, NNLD 2F10.0, 10X, I5 
NVP I5 


VX(1), V¥(1), VZ(1), WX(1), CW(1),DH(1) 6F210.0 


VX(NVP), VY(NVP), VZ(NVP), WX(NVP), 6F10.0 
CW(NVP), DH(NVP) 


The variables on cards 4 and 5 are defined above. In this case, VENTRY 
is only used in determining the force and pressure coefficients and ANG 
is the initial body orientation. The body velocity, wetting factor, 
and increment in depth for each step is defined in cards 7. 


NVP 


VX(I), 
VY¥(I), 
VZ (I) 


WX (I) 


Number of different steps at which entry conditions are 
specified. 


Velocity components in the x, y, z directions of the center 
of gravity (in in./sec.) applied between steps I-1 and I. 


Angular velocity (in degrees/sec.) in the pitch (y-z) plane 
applied between steps I-1 and I. 
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CW(T) Wetting factor applied between the I-1 and I step. If the 
value of this parameter remains constant from step to step 
use the instruction for determining this variable given in 
the CONSTANT orientation section. For the vertical entry 
(VX=VY=0) of pointed bodies an estimate of this parameter for 
each step can be obtained by: 


a. Determining the depth of the entry body, H, below the 
original surface at the start of the step. 


b. Calculating the angle, % between a tangent to the body 
surface and the body axis, z'=H. 


c. Substituting into equation (1) to determine Cy where 


f 


+ 90 - ANG 


ad. For blunt bodies, (nose length/diameter < .75), increase 
this angle by 7% on ogives and decrease it by the same amount 
on cusps. 


DH (IT) Increase in depth (in inches) of the center of gravity 
between steps I-l and I. See instructions in the CONSTANT 
orientation entry section. 


The entry velocity at step I is taken to be the average of that at 
steps I-l and I. It is only necessary to specify data cards for the 
first few steps in which the above parameters change. For steps larger 
than NVP the parameter values at step NVP are used. 


The final set of data cards is used to define the grid on the 
surface of the entry body. The three available options for 
constructing a grid on the body surface are STANDARD, OGIVE, AND LIST. 
These can be used singularly, in combination with one another and can 
be called in arbitrary sequence. The only restriction is that the 
lowest point on the body should occur on that part of the grid 
constructed by the first option called. To indicate the desired 
options, the following input cards are required: 


Card No. Variable Format 
8 N I5 
9.1 option 1 3A4 
9.N option N 3A4 


Here N is the number of options to be used. The recommended options 
are STANDARD and LIST. Description of the OGIVE option will not be 
given in these instructions. For information on this option, refer to 
Reference 19. 
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The grid representing the surface of the entry body should cover 
only the nose of the model and not the afterbody. In all cases the 
pressures on the afterbody are small. Furthermore, on bodies with 
sharp shoulders such as a disk cylinder, the flow separates at the edge 
of the model face. If the afterbody is gridded, the flow is required 
not to separate since the invicid boundary conditions are enforced at 
the centroid of each element. This is physically unrealistic and hence 
neglecting the afterbody is appropriate. 


A description of the three available options follows. Under no 
circumstances should right angles be modeled directly. If the body 
under consideration has such a surface discontinuity, it should be 
modeled with a 89.9 or 90.1 degree angle. 


STANDARD 


This option is applicable to axisymmetric bodies or axisymmetric 
portions of arbitrary bodies. The user specifies rings along which 
nodes are located. Adjacent nodes are combined to form elements. A 
typical grid for a flat, circular plate is shown in Figure G-4. The 
required input is: 


Card No. Variables Format 
10 NROWS, IANG, ISUP 315 
11.1 R(1), Z(1) 2F10.0, I5 
11.NROWS R(NROWS), Z(NROWS), IW(NROWS) 2F10.0, I5 
NROWS Number of grid rings. 
IANG If IANG=0, only half of the face is gridded as shown. 


If IANG=1, the complete face is gridded. 


ISUP If ISUP=1, the stagnation element (element number 1) is 
removed. This option is used for running pointed objects. 
For such bodies, R(1) should be very small (i.e., D/1000) but 
must be finite. 
If ISUP=0 this element is included. 


R(I) Radius of ring I in body fixed coordinates (x',y',z', in 
inches). 

Z (1) z2' coordinate of ring (in inches). 

IW(TI) Number of elements in the area between rings I and I-1l. 


Delete this variable on card 1. If IW=0, elements are 
automatically selected so that they are approximately square. 
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FIG. 6-4 GRID OF A CIRCULAR PLATE 
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WATER 
SURFACE 


LIST 


This option requires that the user input the list of nodes and 
elements to be used in the run and hence is applicable to arbitrary 
bodies. The nodes can be read in any order, however, they are numbered 
sequentially for internal use in the code. Each element is constructed 
using nodes from the input list. The identification numbers of nodes 
defining the four corners of each element must be read in a clockwise 
order with respect to an observer on the outer surface of the element. 
The required input cards are: 


Card No. Variable Format 
10 NP, NE 215 
cies eae x'(1), y'(1), 2'(21) 3F10.0 
11.NP x' (NP), y'(NP), z' (NP) 3F10.0 
12.1 IN(1,1), IN(2,1), IN(3,1),IN(4,1) 4I5 
12.NE IN(1,NE), IN(2,NE), IN(3,NE), IN(4,NE) 415 

NP Number of node points to be read in. 

NE Number of elements to be read. 

x'(I), Location of the Ith node in body fixed coordinates (x', y', 

vr CL) z', in inches). 

z' (TI) 

IN(1,I), Identification number of the nodes defining the four corners 

IN(2,I), of the element. 

IN(3,I), 

IN(4,T) 
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PART 2 


SPECIAL INSTRUCTIONS BEFORE GENERATING FEM 


ENTRY includes three options to define the grid on the surface of 
the entry body: STANDARD, OGIVE, AND LIST. The first two options are 
associated with the automatic grid generation scheme built into the 
ENTRY code and are not suitable for WEST application. The LIST option 
requires the user to input the list of nodes and elements to be used in 
the runstream, and hence is applicable to arbitrary bodies and linkage 
with the PATRAN pre-processor. The translator PATENTR was developed to 
translate the PATRAN neutral file (a PATRAN output ASCII file) defining 
the finite element mesh of the entry body into the water entry model, 
with the same mesh but in LIST option format suitable for use in ENTRY. 
ENTPRES is programmed to convert the ENTRY-calculated element pressures 
and their arrival times into time-varying load specifications in ABAQUS 
format for aynamic analysis and in YADAP format for time-history plots 
on a personal computer. In order to make both translations successful, 
however, several restrictions imposed by ENTRY need to be followed: 


(1) In dividing the surface of an entry body into rings and 
subdividing rings into elements, the node and element identification 
numbers (ID's) must be numbered sequentially, from tip to tail and from 
left to right. If neither begins with the number one (1), a number 
offset must be input for both node and element ID at PATENTR run time. 
The user will be prompted for them at that time. 


(2) Only four-noded quadrilateral (QUAD) elements are permitted in 
ENTRY. The element connectivities which are generally defined in 
counter-clockwise convention in FEA, will be reversed by PATENTR for 
ENTRY. The element pressures calculated by ENTRY are always positive, 
and are treated as positive only when applied in the outward normal 
direction. This direction is determined by the element connectivity 
defined in PATRAN following the right-hand rule. Based on this 
knowledge, the user must decide whether positive or negative pressures 
are to apply to the shell elements in the ABAQUS analysis. The user 
will be prompted for the sign of pressures to be output in an ENTRY run 
when the element pressure-time history output option is activated. 


(3) Define the axis of the entry body as the z-axis. If a FEM is 
already built with x- or y-axis defined as the body-axis, PATENTR can 
convert it into a z~axis model for ENTRY. This is done simply by 
inputting an 'x' or a 'y' when prompted for the body-axis coordinate in 
the FEM at PATENTR run time. Note that this conversion is necessary 
only when running ENTRY. Since pressures are independent of the 
coordinate system chosen, a user can still keep his own body-axis when 
running ABAQUS downstream. 


(4) Although it is not necessary to define the tip of the entry body 
as the origin of the body-axis, it is recommended to do so. This 
action will permit determination of the pressure arrival time for each 
ring of elements to be timed from the tip. This will subsequently 
facilitate the chosing of time-steps in the ABAQUS FEA. Move the 
origin to the tip by inputting a z-axis shift (positive if the tip of 
the ENTRY model is above the origin of the PATRAN FEM) when prompted at 
PATENTR run time. 


(5) It is recommended in the generic use of ENTRY that the increment 
of depth for water entry be chosen such that the average element is 
submerged in two steps. To couple with ABAQUS, however, it is 
necessary to choose the step size precisely so that each ring of 
elements of the entry body will be submerged in exactly one step. This 
is because a pressure cannot be applied to only half an element. Since 
efficient FEM often results in non-uniform mesh density, this situation 
leads to chosing variable step sizes in ENTRY to properly submexge each 
ring of elements in a single step. The times required to submerge the 
variable step-sized rings are computed in ENTRY for given body 
orientation and entry conditions and can be used to synchronize with 
the pressure-arrival time for each ring of elements in an ABAQUS 
analysis. 


(6) The only symmetric condition the ENTRY code can take advantage of 
is half-symmetric, with the y-z plane chosen as the plane of symmetry, 
for either vertical or oblique entry. However, an experienced user can 
take better advantage of symmetry conditions in a subsequent ABAQUS 
analysis by using a quarter- or eighth-symmetric FEM. This becomes 
automatic with the element grouping technique implemented in ENTPRES 
for the case of axi-symmetric bodies entering the water vertically. 


(7) In the case of vertical entry with a zero angle of attack for an 
axi-symmetric body surface, the resulting axi-symmetric pressure 
distribution varies only in the axial direction as the body is 
submerged. Elements in each ring at any given depth will experience 
the same pressure. To facilitate the element pressure-time history 
input for ABAQUS, an element grouping technique is implemented in 
ENTPRES. The user is prompted for the number of elements in a ring 
that have the same pressure at ENTPRES run time. Therefore, it is 
recommended to form a FEM with the same number of elements in each ring 
of the entry body. This grouping limitation excludes the stagnation 
elements at the tip of the ENTRY model. The user will also be prompted 
for the number of such elements included in the model. This exclusion 
increases the flexibility in FEM by allowing the user to use triangular 
(TRI} shell elements at the tip of the body while using QUAD elements 
formed by two TRI shell elements for ENTRY calculation. The near- 
triangle shaped QUAD elements used in ENTRY must therefore be divided 
into two TRI elements for ABAQUS. The grouping technique implemented 
in ENTPRES allows this mixed use of QUAD and TRI elements in ENTRY and 
ABACUS, respectively. 
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PART 3 
COMPUTER RUN STREAM EXAMPLE 


PATENTR, ENTPRES, and options in ENTRY are prompted at run time. 
Options to output pressure-time history files in PATRAN neutral file 
format, options for ENTPRES' conversion to ABAQUS load format for 
analysis, and options to produce YADAP format for time-history plots 
are programmed in an interactive manner. Input/output filenames and 
control parameters are also prompted at run time. Table G-l isa 
sample listing of a VAX 11/785 computer run stream that begins after a 
PATRAN session to generate a FEM of a blunt, 90-degree cone impacting 
water vertically at 100 ft./sec.. Some questions that are annotated 
require explanations: 


(1) CONE90.IN is the file that contains two sets of control parameters 
for water entry. The first set specifies the basic program options for 
ENTRY, such as constant or variable body entry orientation, option to 
print step by step flow field and element information, and the entry 
body's symmetry condition. The second set specifies the water entry 
conditions such as entry velocity, angle of attack, increment of depth, 
and wetting factor. This file is then appended by the PATENTR- 
converted LIST grid file of the cone (CONE.LST) to form the complete 
input data to run ENTRY. The complete iisting of input data for ENTRY 
is given in Appendix A. 


(2) There are two options to input the LIST grid file for ENTRY. The 
first option is appending the LIST grid file to the ENTRY run control 
data file to form a single file for the ENTRY run, as described in note 
(1) above. The second option is to treat the LIST grid file as a 
separate external file during the ENTRY run. Since the first option 
was chosen for this example, the answer to the prompted question is no. 
If the second option is taken, the answer is yes, and the user will be 
prompted for the name of the LIST grid file. 


(3) As shown in Figure 3 of the main text, the inclined surface of the 
cone is divided equally into 10 rings, with 8 QUAD elements each. The 
flat tip is modeled with four near triangle-shaped QUAD elements, as 
shown in Figure 4. 


(4) CONE90.INP is an ASCII file output from running the PATABA 
interface. This file contains all but the loading data required for an 
ABAQUS analysis. It includes the FEM definition of the cone, analysis 
type selection, material and section property cards for elements and 
boundary condition specifications. This file is appended by ENTPRES 
converted time-varying loads (PRES.ABA) for dynamic analysis with 
ABAQUS. A complete listing of input data for an ABAQUS execution is 
given in Appendix Cc. 
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Table G-1. List of VAX 11/785 sample computer run stream. 


$ RUN PATENTR 
ENTER INPUT FILENAME--PATRAN NEUTRAL FILE W/ EXTENSION 
CONE90 .NEU 
ENTER OUTPUT FILENAME--ENTRY "LIST" GRID OPTION FILE 
CONE90.LST 
IS NODE OR ELEM NUMBER SHIFT REQUIRED FOR ENTRY? Y/N 
IF YES, ENTER NODSFT AND NELSFT IN 2I5 FORMAT. 
NOTE: NODE & ELEM ID MUST BEGIN WITH 1 & IN SEQUENCE 
N 
ENTER BODY-AXIS CHOSEN IN FINITE ELEM MODEL, X,Y OR 2 
NOTE: NON-Z SYSTEM WILL BE CONVERTED TO Z-X-Y SYSTEM 
WITH Z AS BODY-AXIS & Y-Z AS PLANE OF SYMMETRY 
Z 
IS AXIAL OFFSET REQUIRED FOR ENTRY MODEL? Y/N. IF YES, 
ENTER SHIFT (POSITIVE IF TIP IS ABOVE ORIGIN OF FEM) 
NOTE: THIS ALLOWS WETTING TO BE TIMED FROM THE TIP. 
N 
END OF INPUT FILE ENCOUNTERED--NORMAL TERMINATION 
FORTRAN STOP 
$ 


$ APPEND CONE90.LST CONE90.IN (1) 
$ RUN ENTRY 
ENTER INPUT FILENAME INCLUDING EXTENSION 
CONE90.IN 
ENTER OUTPUT FILENAME INCLUDING EXTENSION 
CONE90.OUT 
ENTER OUTPUT FILENAME FOR ELEMENT PRESSURE-TIME 
HISTORIES IF WANTED, TYPE NO IF NOT WANTED 
CONE90.PRE 
ENTER FILENAME FOR TOTAL FORCE-TIME HISTORIES, DRAG & 
NORMAL FORCES & MOMENTS @ CG IF WANTED, TYPE NO IF NOT 
CONE90.F 
IS LIST INPUT DEFINING GRID FROM AN EXTERNAL FILE? Y/N 
N (2) 
FOLLOWING ENTRIES COMFORM TO PATRAN NEUTRAL FILE FORMAT 
FOR DISTRIBUTED LOADS: 
(1) TITLE CARD (DATA PACKET TYPE 25) 
(2) SUMMARY DATA (DATA PACKET TYPE 26) 
(3) DISTRIBUTED LOADS (DATA PACKET TYPE 6) 
ENTER TITLE OF PRESSURE-FILE (LIMIT 80 CHARS) 
PRESSURE-TIME HISTORIES OF 100 FPS VERTICAL ENTRY BLUNT 90-D CONE 
ENTER 6 LOAD COMPONENT FLAGS (ICOMP=0/1) AT 611 FORMAT 
EX: FOR PRESSURE LOADING ON QUAD ELEM, ENTER 001000 
001000 
ENTER 8 ELEM NODE FLAGS (NODE=0 OR 1) AT 811 FORMAT. 
1 MEANS THAT NODE IS ON THE LOADED EDGE OR SURFACE. 
EX: FOR A PRESSURE-LOADED QUAD ELEM, ENTER 11110000 
11110000 
ENTER 1ST LOAD SET ID (IV=1) & ELEM FACE NO. (NFE=1-6) 
AT 215 FORMAT. NFE=0 FOR A PRESSURE-LOADED QUAD ELEM 
1 e) 
ARE PRESSURES APPLIED TO PATRAN MODEL +(P) OR -(N)? P/N 
NOTE: PRESSURE IS + IN POSITIVE NORMAL FOR QUAD ELEMS 
N 
FORTRAN STOP 
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Table G-1. (continued) 


$ RUN ENTPRES 
ENTER INPUT FILE OF ELEM PRESSURE HISTORY, EX: CONE.PRE 
CONE90.PRE 
ENTER OUTPUT FILE FOR ABAQUS FE ANALYSIS, EX: PRES.ABA 
PRES. ABA 
DO YOU WANT TAPE9.DAT TIME-HISTORY PLOT FILE? Y/N 
NOTE: TAPE9.DAT FORMAT IS 12E11.4, WITH ELAPSED TIMES 
IN FIRST COLUMN AND UP TO 11 GROUP PRESSURES 
THEREAFTER. IF MORE THAN 11 GROUPS OF OF DATA 
ARE ENCOUNTERED, THEY WILL BE SPREADED OVER TO 
FORO10, FORO11, FORO12.DAT...IN THE SAME FORMAT 
Y 
ENTER TITLE FOR TAPE9.DAT PLOT FILE (80 CHARS MAX) 
PRESSURE-TIME HISTORIES OF 100 FPS VERTICAL ENTRY BLUNT 90-D CONE 
ENTER NOS. OF STEPS (NSTP,I5) AND NOS. OF ELEMS IN A 
GROUP (NEL,I5) WITH THE SAME PRESSURE IN INPUT FILE 
NOTE: IF ALL ELEMS HAVE DISTINCT PRESSURES, SET NEL=1 
it 8 (3) 
ARE STAGNATION ELEMS INCLUDED IN WATER-ENTRY MODEL? Y/N 
Y 
(3) 
IF YES, ENTER NUMBER OF STAGNATION ELEMS (NSTG,I2) 
NOTE: STAGNATION ELEMS ARE TREATED AS GROUP NO. 1 
SUMMARY OF PROCESS: 
MAX NUMBER OF PRESSURIZED ELEMS WITHIN STEPS 
NUMBER OF ELEM GROUPS (WITH EQUAL ELEM PRESSURE) 
NUMBER OF ELEMS PER GROUP (EXCL STGN ELEM GROUP) 
NUMBER OF ELEMS IN STGN ELEM GROUP (GROUP NO. 1) 
NUMBER OF TIME-HISTORY PLOT FILES GENERATED 
FORTRAN STOP 


oid wou 
[ee] 


$ 
$ SET DEF [RCSHAW] 
$ APPEND [.ENTPRES]PRES.ABA CONE90.INP (4) 
$ @ABAQUS 
IDENTIFIER : 
INPUT FILE NAME (W/O .INP) : 
RESTART READ FILE ? 
OLD ‘FILE OUTPUT FILE’ FILE i : 
BATCH QUEUE YOU WANT (SYS=<CR> SLOW=S) ? : 
GIVE THE AMOUNT OF TIME YOU WANT TO HOLD THE JOB FOR 
THE FORMAT IS HH:MM: ( <RET> TO IGNORE) 
CONE (queue SYSSBATCH, entry 350) started on SYS$BATCH 


Job CONE (queue SYSS$BATCH, entry 350) completed 
$ DIR CONE 
Directory USER2: [RCSHAW] 


CONE .COM;1 CONE.DAT;1 CONE.FIL;1 CONE. LOG; 1 
CONE.STA;1 


Total of 5 files. 
$ 
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